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FOREWORD 


IRE Standards on Symbols and Terminology for Feedback Control Systems 


Although the standards on symbols and terminology were published in the Nov- 
ember, 1955, and January, 1956, issues of the Proceedings, they are reproduced in 
this Transactions to remind control engineers that standards do exist. Standard 
symbols and terminology definitely increase the efficiency of obtaining informa- 
tion from technical material. Unfortunately, all papers in this issue of the 
Transactions do not use the standard symbols, but an attempt will be made in the 
future to have authors use standards as much as possible. 


Considerable debate, argument and compromise was necessary to arrive at 
these standards; each term was carefully chosen. Most of the terminology is com- 
patible with other proposed terminology in the field. However, it is intended 
to be more fundamental, applicable to a general feedback loop as well as for a 
specific form. For example, the loop input signal may be applied anywhere in the 
loop; the loop output signal may be chosen anywhere in the loop. 


It was felt that this was necessary to define adequately the characteristics 
of the teterogenous control loops, network and systems which are encountered by 
the diversified interests of the IRE professional groups. Any criticism of these 
terms or suggestions for modification will be seriously considered by Subcommit- 
tee 26.1 and Committee 26.0, and it will be discussed in the Transactions. 


Meanwhile, careful study and use of these standard symbols and terminology 
by members of the PGAC will simplify the problem of conveying new ideas with a 
minimum of effort. 


Abstracts of Minutes of the Meeting of the Administrative Committee 


The committee felt that PGAC members might desire to know the problems, in- 
terests and scope of the group as a whole. The minutes of the Administrative 
Committee meetings are to be abstracted and included in the Transactions for this 
purpose. Anyone interested in serving on the Administrative Committee or nation- 
al subcommittees should write to the PGAC chairman or the editor of the Transac- 
tions. Assistance on future committee work will be needed. 


Index of PGAC Membership 


The membership of the PGAC has increased steadily since its inception. To 
illustrate its growth and to indicate the membership in various sections, the 
names of all the members are given in this issue. This should be an aid to those 
wishing to form chapters which will strengthen the PGAC and enrich the Transac- 
tions with papers from local meetings. 


It is hoped that more IRE members interested in control work will join the 
PGAC and identify themselves professionally in this automatic control index. 
The index will be published periodically in the future. 


THEORY, PRACTICE OR INSTRUCTION? 


The first four papers in this issue are essentially mathematical in nature. 
Although examples of applying the theories are given, there are no actual mea- 
surements or applications to physical equipment correlating the theory with 
practice. In general, the ideas should be contributions to the field of auto- 
matic control, and it is hoped that they may be useful in solving practical 
problems or in stimulating other unique methods of analysis or synthesis. Papers 
on the application of these and other theories to control systems will be partic- 
ularly welcome. 


It is interesting to note that responses of the nonlinear systems are given 
exactly, without recourse to the popular describing function or phase plane 
methods. In fact, one nonlinear system is synthesized in such a manner that it 
becomes a relatively simple system which is easy to analyze exactly. Perhaps 
this idea could be applied to a larger class of systems. As pointed out in one 
paper, exact solutions can be used to check the current approximation techniques 
and others which may be developed. 


The last two papers in this issue are tutorial. Relatively short tutorial 
papers are basically difficult to write. From the reader's viewpoint the back- 
ground material may be too extensive or too brief, the subject matter may be too 
sketchy or unbalanced, some pertinent facts may be omitted. All of the field 
cannot be covered, and few concrete examples can be given. The last paper on the 
survey of the sampled-data techniques contains conclusions about methods of anale- 
ysis which undoubtedly would be modified by some readers more familiar with other 
techniques. Comments on these tutorial papers and suggestions for other topics 
will be welcome. Actually, reactions to the first PGAC Transactions were favora- 
ble, but few. More problems, discussions, notes, opinions and comments on 
improving future issues would be desirable. In particular, it would be interest- 
ing to have opinions as to whether papers on theory, practice or instruction 
should be predominant. 


Of course the Transactions are issued for those interested in automatic 
control and in advancing knowledge in the field. This can be done most effec- 
tively by using the Transactions as a medium to interchange ideas, stimulate 
discussions and define the issues which confront the automatic control engineer. 
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THE ISSUE IN BRIEF 


FINAL VALUE CONTROLLER SYNTHESIS ...esecccvccccceeee Ms Ve Mathews and Ce. We Steeg 


A method is presented for synthesizing a type of control system designated as a 
final time controller. This device is a feecback control system that is designed 
to achieve a desired response at one time only, the response at earlier times 
being arbitrary within physical limits. In addition to a dynamic element which 
has a response that is to be controlled, the control loop consists of a feedback 
component and a controller, both of which have characteristics derived by a syn- 
thesis procedure. 


The synthesis procedure presented results in a time varying, nonlinear sys- 
tem sufficiently simple so that optimization with respect to Gaussian random 
disturbances and to initial transients may be made easily. A time nonlinear 
optimum system is achieved, where the nonlinearity is saturation, the principal 
limitation for most control systems. 


Application of the procedure reduces the analysis of the closed-loop system 
performance to the investigation of a single first-order differential equation 
involving the impulse response of the controlled element. 


A POSITIONING SERVOMECHANISM WITH A FINITE TIME DELAY 
AND A SIGNAL LIMITER e@eeeeevseoeveeceeeeeveeev een eeoeneeeeeoeeeeeneoonoeevneoeaeoeeveasee De te Evans 


An idealization of a servomechanism which is used in a digital positioning cir- 
cuit is analyzed. The loop consists of an ideal integrator, a finite time delay 
and a nonlinear signal limiting element. The nonlinear element in the physical 
system contains logical circuitry which converts the digital information into a 
continuously variable error. 


An exact solution using contour integration is obtained for the idealized 
system to provide a standard with which solutions obtained by approximation 
techniques may be compared. In addition, a simple approximation solution and 
an error term are obtained. 


Finally, curves are given which express the time necessary to zero in on ~ 
the final position from a given initial position as a function of the loop param- 
eters. These curves indicate that a minimum, noncritical settling time may be 


- 


obtained by proper selection of loop components. 


THE STABILIZATION OF NONLINEAR SERVOMECHANISMS ENCOUNTERED 
IN ANT=NNA INSTRUMENTATION CCT HOSTEOCHEOHRSHHTHSOSSCHOSHHHRHSHSHSEHOHHHHHSCHHHEHHHOHHE Je Bacon 


This paper concerns the problem of compensating the loop-gain of instrument 


‘servos to provide uniform transient response over a wide dynamic operating range. 


Particular consideration is given to applications in the realm of antenna 
measurements. 


A distinction is made between situations where the nonlinearity results from 
characteristics of the follow-up device and where it results from the loop-gain 
being functionally related to an external variable. In the former case the loop- 
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gain is made self-linearizing; in the latter, the gain is effectively equalized 
with an auxiliary logarithmic servo. This necessarily requires a suitable slave 
potentiometer in nonlinear loop, actuated by the logarithmic servo. 


Certain advantages accure from using a ladder attenuator as the slave com- 
ponent. These are described. It is shown how this combination logarithmic servo 
and ladder slave attenuator can make a loop-gain invariable which otherwise would 
be functionally related to a variable E by the equation G = EEN, The idea is 
extended to the generation of polynominal terms of the same form as G. 


ON THE DESIGN OF A=C NETWORKS FOR SERVO COMPENSATION «ccccovcceeeee He Levenstein 


An analytical method for analysis and synthesis of networks for AC servo compen= 
sation is presented. The response of a linear network to a modulated suppressed-= 
carrier excitation is formulated in terms of in-phase and quadrature carrier com- 
ponents. The relationship between the modulation on these components and the 
exciting modulation is shown to depend upon operators simply related to the orig- 
inal network function. 


The expansion of the data-frequency operators into partial fractions is 
shown to lead to a simple synthesis procedure for deriving the original network 
operator from the in-phase or quadrature operators. As an example, the process 
is applied to the derivation of a representative network for lead compensation 
of an AC servo responsive to the in-phase component of error signal. The use of 
RC networks in AC servo compensation is shown to be limited to derivative types 
of equalization. 


FUNDAMENTAL EQUATIONS FOR THE APPLICATION OF STATISTICAL 
TECHNIQUES TO FEEDBACK CONTROL SYSTEMS .ccoccccccccccececceceee Ge As Biernson 


Basic equations necessary for applying statistical techniques to the design of 
control systems are presented clearly and concisely to make them more understand- 
able and useful to the engineer. 


A discussion of the difference between noise and signal inputs and the need 
for ‘statistical techniques is followed by Gaussian expressions for noise ampli- 
tude density, the probability of noise saturation of system components and prac- 
tical considerations in determining proper components saturation levels. After 
describing autocorrelation functions of input and output, the corresponding 
cross-correlation functions and their relation to rms values are derived, and the 
relation between statistical inputs to a system and the resulting rms output is 
developed by transient and by transform methods. 


Appendices show graphical interpretations of convolution and correlation 
integrals, the relations between spectral density and time function transforms 
and a method for using an analog computer to generate mean-square values of sys- 
tem outputs for statistical inputs. 


A SURVEY OF TECHNIQUES FOR THE ANALYSIS OF SAMPLED=DATA 
CONTROL SYSTEMS welpieie 016 6W4 «els 0's u'e'We slainca he's iain ciel a taeri ts Murphy and Rel. Ormsby 


The present use in control systems of pulsed-data links, track-while-scan radar 
and digital computers has stimulated interest in the analysis of sampled-data 
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control systems. The effort now being expended to develop techniques for 
analyzing such systems has resulted in several methods of analysis. 


This paper presents a discussion of three principal methods of analysis 
currently used: the impulse response, frequency response and Z transformation 
of the system transfer functions. The methods are explained and applied, and 
the advantages and limitation of each are discussed. 
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FINAL VALUE CONTROLLER SYNTHESIS 


M. V. Mathews 
Bell Telephone Laboratories, Inc. 
Murray Hill, New Jersey 
and 
C. W. Steeg 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 


Summary -- A method is presented for synthesizing a type of control system desig- 
nated as a final time controller. This device is a feedback control system that 
is designed to achieve a desired response at one time only, the response at ear- 
lier times being arbitrary within physical limits. In addition to the dynamic 
element, which has a response that is to be controlled, the control loop consists 
of a feedback element and a controller, both of which have characteristics 
derived by the synthesis procedure. The final time controller developed is a 
time varying, nonlinear system and the procedure is one of the few that are 
applicable to the synthesis of nonlinear controllers. Application of the pro- 
cedure reduces analysis of the performance of the closed-loop system to the 
integration of a single, first order differential equation involving the impulse 
response of the controlled element. 


DESCRIPTION OF FINAL VALUE CONTROLLERS 


This paper presents a method for synthesizing a class of feedback control 
systems, which differ from ordinary servomechanisms in that the controlled system 
response or output is to assume a desired value at one particular time only, the 
response at earlier times being arbitrary. Systems of this type are designated 
here as final value controllers. An example of final value control is return-to- 
base operation of an aircraft. Efficient procedure requires that the aircraft 
approach its landing field in such a way as to pass through a ground-control- 
approach (GCA) gate at a particular time but the trajectory flown at previous 
times can vary within wide limits. 


Many of the existing methods for the design of systems whose effectiveness 
depends — on satisfying final conditions utilize ordinary servomechanism tech- 
niques’~’ which attempt to zero an error signal at all times. However, the 
unnecessary restriction of minimizing an error weighted equally for all times 
may unnecessarily complicate the design procedure and reduce etka tol) per- 
formance achievable. 


The synthesis procedure presented here results in a sufficiently simple 
time varying, nonlinear system so that optimization with respect to Gaussian 
random disturbances and to initial transients may be made easily. A true non- 
linear optimum system is achieved, where saturation in the prime mover is the 
nonlinear limitation. This factor is the principal nonlinear limitation for 
most control systems. 


For the purpose of explanation, one particular final value control situation 
is considered and is idealized into the system shown in Fig. 1. For this system 


6 


INPUT CONTROLLER DYNAMIC 
PREDICTOR F (Xip» Xrp, t) ELEMENT 


CONDITIONAL |__| 
RESPONSE 
PREDICTOR 


INFORMATION 

DESCRIBING 

*“STATE” OF 
DYNAMIC ELEMENT 


Fig. 1 - Final value controller. 


the final value problem will be defined, the operation of the system described 
physically and the optimum characteristics pointed out. 


STATEMENT OF THE PROBLEM 


The problem consists of so synthesizing the several elements of the control 
system that the final response of the dynamic element x,(T) will equal some 
desired value xq, the time T being the final time. The desired value xq is 
related statistically to the input x3 (t) with the result that xq can be predicted 
from x;(t); e.g., x(t) might be xq plus noise. The given data for the problem 
consist of the input x;(t) and the dynamic element that is to be controlled. The 
synthesis process gives specifications for the operation of the input predictor, 
the controller and the response predictor through minimization of the final error 


E where 
E = xq - x,(T). (1) 


The dynamic element is assumed to saturate at its input but otherwise to be 
linear. Thus, for the magnitude of x,(t) less than the saturation value xg, the 
response xy(t) is related to x,(t) by a linear differential equation and for 

ko (t)] > Xom the effective input to the linear portion of the dynamic element is 
<xom- This statement constitutes the definition of the types-of dynamic elements 
That are considered. | 


The physical operation of the system can be described as follows: at any 
time t the input predictor examines all past values of the input x; and from 
these values makes the best prediction of the final value xg. This prediction is 
designated x;,(t). At this same instant of time the conditional response predic- 
tor tells whad will be the final value of the response if the control signal from 
t till the final time is some known function gyn. This conditionally predicted 
response is designated xrp(t). If xrp(t) equals Xip(t), then the controller acts 
to make gn the control signal. If xpp(t) equals xjp(t), then the controller 
applies the maximum control signal consistent with saturation in the dynamic | 
element in order to reduce the predicted terminal error. : 


Two properties achieved by the system are: - (1) the prediction of xg from 
x; (t) is completely separated from the closed-lcop control of the dynamic element 


i. 


and (2) the closed-loop=control equation is reduced to a first order differential 
equation no matter what the complexity of the dynamic element. 


The prediction of xg from s(t) and the synthesis of an ideal input predic- 
tor has been treated by Booton.4” He shows that the system achieves a minimum rms 
error when the random component in the input is Gaussian. This paper is con~ 
cerned with the closed-loop=control equation and the trensient performance. 
Rooton's development will not be repeated. However, one characteristic of the 
input predictor is necessary to justify the choice of gn, the equilibrium control 
function, and to. justify the transient performance criteria. If the input pre- 
dictor is ideal x; (t) will consist of xg plus a completely unpredictable noise. 
Consequently, because x;,(t) has an equal probability of increasing or decreas- 
ing, the best function gn is zero, so that the full positive or negative control 
will be available to correct for unexpected changes in x (t). This hueristic 
argument can be made exact when x;(t) equals xg plus a Gaussian random noise. 


Because the variation in Xip(t) consists of an unpredictable noise, there 
can be no predictable transients in Xi p(t) such as terms which vary as t or t*. 
Consequently, the only significant transient response is that to an initial step 
change in xj)(t). 


The principal assumptions and simplifications which have been used in the 
system of Fig. 1 are the following: 


1. All statistical uncertainty is associated with x;(t). The response and 
state of the dynamic element are assumed to be exactiy known so that if 
the control input were gn, the conditionally predicted response would be 
exactly achieved. This assumption amounts to neglecting uncertainty about 
the nature of the dynamic element and noise in the measurement of its 
state. In some cases this uncertainty can be transferred to x;(t). 


2. Only one final error, E, is minimized. Systems which simultaneously mini- 
mize a number of final conditions can also be developed. 


3. The dynamic element is assumed to saturate at its input as previously 
discussed. 


; The action of the closed-loop system which has been described above is 
analyzed in the next section. ; 


CLOSED=LOOP RESPONSE EQUATION 


To develop the differential equation which describes the performance of the 
closed-loop system the operation of the controller and the conditional response 
predictor must be characterized. 

The controller may be specified by a nondifferential function F, 

Xo (t) =F [xip(t), Xpp(t), t | (2) 


which relates x,(t) to x4p(t), x,p(t), and t. The function will be in such form 
that x,(t) never exceeds Ene saturation limit xgme Usually F will have the 
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properties of an off-on relay type of controller as will be shown. In order tu 
determine Xpp(t), the final response may be written in the form 


we 
x,(T) = | ny Efe Gd sg (3) 


- 00 
in which h(t) is the impulse response of the unsaturated dynamic element. Equa~ 
tion (3) applies only for a constant-coefficient dynamic element. However, the 
generalization to a time varying element can easily be made. For simplicity, 
only the constant-coefficient system is discussed. The integral in Eq. (3) 
extends from minus infinity to T. Because the controller is started at some 
finite time Tj, the integral from -oo to Tj implicitly specifies the initial 
conditions at T). 


At time t the conditional response predictor assumes that x, (t) becomes and 
remains a known function g, which in the present example is taken to be zero. 
Thus, because x,(1) is taken as zero for t<T<T, the conditional response 
xp(T) cyt which is also the conditional predicted response may be written: 


t 
xr(T)|¢,¢ = Xp(t) = af h(T-7) x,(T)dT. (4) 
=0O 
The differential equation for the closed-loop system is obtained by 
differentiating Eq. (h) to yield 


rp = h(T-t) x(t) (5) 
at 
and combining this relation with Eq. (2) to yield 
Srp - n(t-t) Flxpp(t), xpp(t), t] = 0 (6) 
dt 


‘which is the desired equation describing the performance of the closed-loop 


system. Equation (6) can be integrated to obtain xp. Since at the final time, 
Xpp(T) = xy(T), the error E may be written 


E = xq - X,p(T) | (7) ) 


and a solution for E may always be obtained from Eq. (6). Usually this error is 
sufficient to evaluate a controller. However, the complete behavior of the sys- 
tem can be computed from xyp(t). Thus, Eq. (6) completely describes the final 
value performance of the system. In addition, it is a simple, first order 
equation. 


Most ‘of the desirable properties of final value controllers arise from the 
first order character of Eq. (6). Some of these properties can be illustrated 
by specifying a particular control function as shown in Fig. 2 and given by the 
equation . 


\ 


- Kee 3 + saturation region 
( t 
Xe(t) = K [ip() - Xp (t)} linear region (8) 
‘i : ( -x ang - saturation region. 
‘cm 


When the gain K becomes infinite an initial transient is removed in a minimum 
time and a minimum rms error due to Gaussian noise is achieved. The transient 
properties will now be developed. 


In the absence of noise and with an ideal input predictor x, p(t) is constant 
and equal to xq. yes the transient behavior is determined solely by the time 
necessary for (t) to change from its initial value to some specified constant. 
This time is ae cates most directly by Eq. (5). If the predicted terminal error 
is nonzero x,(t) will be +x,, and 


Crp = h(T=t) (xem). (9) 


h(t) > O« (10) 


for all values of t then Eq. (9) can be interpreted as showing that x_ (t) 
increases or decreases monotonically at the maximum rate as limited EP the 
saturation level x,,, and the impulse response of the dynamic element h(t). 

Thus the initial error is reduced to zero for a minimum value of T as compared 
with any other possible controller. 


If K is finite and x,,, sufficiently large that the system operates linearly 
then in the transient case Bas (6) becomes 


dxrp + Kh(T-t) (t) = Kh(T-t) x 11 
which can be integrated explicitly and the error evaluated at the final time to 
give, T 
E = Xq - Xrp(T) = [xg - Xpp(o) | exp [-« J h(t) dt | (12) 


in ee the quantity [x, - Xpp(o)] is the final error if no control has been 
appslLede . 


*If h(t) does not satisfy condition (10) then the same optimum behavior can be 
obtained by a controller whose gain changes sign each time h(T~t) changes sign. 


F (Xip, Xrp,t) 


Fig. 2 = Nonlinear controller. 
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Equation (12) shows that if condition (10) is again satisfied the terminal 
error will decrease monotonically as the gain K is increased and as the terminal 


time T is increased. Thus the controller can be said to be stable for any value 
of gain. 


This section has pointed out that the essential performance of the control 
system can be calculated by integrating one first order differential equation, 
Eq- (6). In addition, the optimum transient performance and stability of the 
system was indicated by an examination of the differential equation with one 
particular control: function. Thus, the essential properties of the final time 
controller have been put in. evidence. There remains to discuss the form of the 
conditional response predictor and present an illustrative example. 


SYNTHESIS OF “THE CONDITIONAL RESPONSE PREDICTOR 


The function of the conditional-response predictor has been specified in the 
preceding sections. At time t the predictor must compute the response of the 
dynamic element at the terminal time T under the condition that the control sig- 
nal xe(t) is zero. Because the control signal is assumed to be zero during the 
interval from t to T, the response prediction is simply an evaluation of a homog- 
eneous solution of the dynamic element differential equation.* This homogeneous 
solution must satisfy at time t a set of initial conditions which are by defini- 
tion the state of the dynamic element. A set of state variables are defined as 
those that are sufficient to specify all the arbitrary constants in the complete 
homogeneous solution. For example, if the dynamic element differential equation 
is of the form 


an xf") (4) + an] x(n-1) (+) + coe + Ag Xp(t) = Xe(t) (13) 


then the response xy(t) and its n-l derivatives x1) (+4), ate xfn-1) (+4) make up a 
complete set of state variables. For this example, the conditional response 
prediction is most simply written in terms of the homogeneous solutions 

f£3(t), i = 0, 1, ..., n-l which have the initial values 


ef (o) = (1a ej. (14) 
(0,1 # 3 


With these functions the homogeneous solution having the state variables as its 
initial condition at time t is simply 


n-1 . 
F(u) = => xr(4)(t) ty(u = t) . (15) 
1=0 
In Eq. (15) u has been introduced as a new independent time variable to avoid 
confusion with t. That F(u) satisfies the required initial conditions can be 
established by taking its j h derivative and setting u = t which as a result of 


Eq. (14) yields 


het G4) 9 rah (4) ote 
PD (t) = = x) (+) tf," (0) = xp/(t) . (16) 


*In the case of a nonzero gp function the prediction would be the sum of the 
homogeneous solution and a known particular solution. 
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The conditional response at the terminal time is simply F(T) and thus 
n=1 
i 
F(T) = xpp(t) = =, x2) (t) 24 (7-4) . (17) 


Equation (17) can be interpreted physically to synthesize the predictor. 

The response Xp (+) and its derivatives form the inputs to the predictor and the 
functions f;(T-t) are generated within the predictor. The entire predictor is 
realized as shown in Fig. 3. If the dynamic element differential equation is 
more complicated than Eq. (13) then an expression as simple as Eq. (16) may be 
obtainable for the conditional response predictor. However, the general proce= 
dure of finding a homogeneous solution and a complete set of initial conditions 
can always be used to construct the predictor. 


The discussion of the example gives the false impression that cnt) and its 
derivatives are usually measured as state variables in a physical element. This 
procedure is seldom desirable since the higher derivative measurements will be 


Fo (T-t)  # (T=t) fy -; (T-t) 


Fig. 3 - Conditional response predictor. 


affected by noise. Instead, for almost all cases a set of easily measurable 
state variables can be determined as the Sper ey. stored in various forms in the 
ed Sata element. 


ANALYSIS OF AN EXAMPLE , 


A realistic. yet elementary example with which to illustrate the final value 
control theory is difficult to obtain. The theory was developed too recently for 
many systems to have been examined and most of the applications deal with compli- 
cated controllers which satisfy more than cne terminal condition. As a conse- 
quence a simple hypothetical system will be presented, the conditional response 
predictor equations developed and the transient response computed. 


The control system is shown in Fig. and consists of the usual elements: 
input predictor, controller, dynamic element and conditonal response predictor. 
The dynamic element saturates at a control-signal magnitude Xeme The unsaturated 

element consists of a simple time constant followed by an integration and is 
des¢ribed by ‘the equation 


\ 


Txh (t) + xp (t) = xQ(t). panmectice' 27, 


a 


1 
1“ 


CONTROLLER 
= 


INPUT 
PREDICTOR 


CONDITIONAL RESPONSE 
PREDICTOR 


Fig. h - Simple final value controller. 


The controller is assumed to be the ideal controller shown in Fig. 2 with an 
infinite K. The conditional response predictor equations are developed below. 
We will assume the input predictor to be ideal so that to evaluate the transient 
performance in the absence of noise X4 p(t) = Xge The impulse response of the 
dynamic element may easily be shown to be 


h(t) =1-07% a9) 
and the homogeneous response of the dynamic element to be 

CE yiCek Gage (20) 
where A and B are arbitrary constants of integration. 


To determine the conditional response predictor equations the state variable 
must be specified first. These variables are simply x,(t) and its derivative 
x(t). The homogeneous response of the dynamic element with initial condition 
xi(t), x,(t) can either be computed from a direct analysis of the dynamic element 
or from Eq. (20). With either procedure the response f,(t) corresponding to 
x,(0) = 1, x8(0) = 0 is 


fo(t) = 1 (21) 


and the response f}(t) corresponding to x,(0) = 0; x1(0) = 1 is 


a | 
f,(t) =t(l-e%). | (22) 


Thus, the conditional predicted response at time T corresponding to x,(t), xp (t) 
at time t is 


y T= 
xp(T)| 0,4 * Xrp(t) = x(t) + xE(t) EE- a oe (23) 


The conditional response predictor shown in Fig. ) is an instrumentation of 


Eq. (23). Thus, the synthesis of the control system is complete. 
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The terminal error can be evaluated as xq = xrp(T) by integrating Eq. (24) from 
0 to T to yield 


alin 


E ® Xq - X(T) = Xq - (T + Te 2410 bi rae Se (25) 
Equation (25) specifies the terminal error E for values of T less than some 
critical value T, where the error becomes zeroe For T>T, the final error is 
always zero. 


A plot of normalized final error E/xg is shown in Fig. 5 as a function of 
normalized time T/t and for several values of saturation xg,- In Fig. 6 the 
final value control system performance is compared with that of a well—adjusted 
linear servo using the same dynamic element. The block diagram of the servo is 
also sketched in Fig. 6. The ratio of Xq to X47 has been so selected that the 
maximum value of xc in the linear servo just equals x,,- As demonstrated in 
Fig. 6, the final value controller not only achieves a zero error for a smaller 


The transient response can be evaluated by assuming the conditions 
xy(0) = x4(0) = 0, and x;,(t) = xq, The final error can then be computed as a 
function of the final time T. The error will decrease to zero as T is increased. 
In the region of nonzero terminal error xq - Xpp(t) will also be nonzero, 
O<t<T, and the controller will apply the maximum input xg, to the dynamic 
element during the entire trajectory. Thus, from Eq. (9) 
T=t 
Press oF TEA) S (2h) 
Pata cm 


value of T than required by the linear servo but also the error remains zero for 
T>T,., while the linear system has a substantial overshoot. 


CONCLUSION 


A procedure has been developed for the synthesis of certain time varying, 
nonlinear systems. The theory tends to show the advantages of these systems over 
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Fig. 5 - Final value controller performance. 
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Fig. 6 - Comparison of linear servo and final value controller. 


systems designed by previously published methods. In some cases this synthesis 
procedure results in the best possible performance obtainable with a particular 
desired response and with a given dynamic element. The procedure is an example 
of a synthesis involving nonlinear differential equations and is simpler than 
most analytical methods for such equations. This simplicity is achieved because 
the synthesis procedure is chosen to result in an analytically tractable system 
(specifically to result in a first order system). Consequently, the requirement 
of tractability has been added to the usual synthesis requirement of good system 
performance. Tractability may be the factor that has limited the use of nonlin- 
ear synthesis schemes. 


Several direct extensions of the procedure which were not discussed are 
apparent. The final value controllers were designed to satisfy only one terminal 
condition, namely, that the response of the dynamic element at the terminal time 
equal some desired value. Systems that satisfy additional requirements involving 
the response and its derivatives at one or more terminal times can be developed. 
The assumption was made that the equilibrium value of the control input to the 
dynamic element is zero. Other controllers can be derived in which the equilib- 
rium value can be any function of time. 


The synthesis procedure is applicable to dynamic elements which are linear 
but saturable. While these elements undoubtedly form a very important class, the 
synthesis procedure probably can be extended to include other types of dynamic 


nonlinear elements. 
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A POSITIONING SERVOMECHANISM 
WITH A FINITE TIME DELAY AND A SIGNAL LIMITER 


David H. Evans 
Bell Telephone Laboratories, Inc. 


Summary -- An idealization of a servomechanism which is used in a digital posi- 
tioning circuit is analyzed. The loop consists of an ideal integrator, a finite 
time delay and a nonlinear signal limiting element. An exact solution is 
derived. Also, a simple approximate solution and an error term are obtained. 
Finally, curves are given which express the time necessary to zero in on the 
final position from a given initial position as a function of the loop 
parameters. 


INTRODUCTION 


There are two salient features about the model to be analyzed -= the nonlin- 
ear signal limiting element and the finite time delay. Although the physical 
system which led to this analysis will not be considered we will indicate the 
sources of these features. 


The nonlinear element derives from the fact that the physical system is a 
digital positioning servomechanism. The nonlinear element in the physical system 
contains logical circuitry which converts the digital (and transitional digital) 
information into a contimously variable error. The model actually has the non- 
linear element at a different place ‘in the loop than is indicated by the above; 
this was for ease of analysis and is immaterial to the applicability of the 
analysis. 


The finite time delay arose from the high speed at which the physical posi- 
tioning servomechanism operates. The positioning speed relative to the finite 
speed of propagation of signals through the system made the delay due to signal 
propagation noninfinitesimal. In the model the delays are lumped. 


THE LOOP EQUATIONS FOR THE MODEL 


In the model shown in Fig. 1. the feedback loop has unity gain. The box 
labeled A contributes a delay of A sec; that is, if the input is v(t+A), the 
output is v(t). f is a nonlinear device such that if v(t) is the input f[v(t)] 
is the output. A is an amplifier. 6,;(t) and 0,(t) are the input and output, 
respectively. 


A = OELAYA SECS 
f = NON-LINEAR ELEMENT 


For a step function input of amplitude M at time t = 0, and 05 (0) S10; eee, 
an initial error M, 


t (1) 
u(t +A) =~ f ffu(x)] K(t - x) dx, t>0 
ie) 


u(t) = M O0<t<A 
where 
u(t) = Q, (t) - 0 (t) = error 
K(t) = 8-function response of amplifier. 
The particular nonlinear f which we use is shown in Fig. 2. Analytically, 
y, u(t).>> 
f[u(t)] = ~ u(t), |u(t)| an 

-y, u(t) < =r 


t [uct] 


u(t) = ERROR 
t[u(t)] = ERROR SIGNAL 


Pirsuc 


The particular amplifier we wish to consider is one with a 6-function 
response 


K(t) = ape P*, 


K(t) can be approximated very closely for this problem by an ideal integrator. 
Hence, take 


K(t) = AJB. | (2) 


With the above assumptions a qualitative picture of the response is shown 
in Fig. 3. The diagram illustrates why we can use an integrator to approximate 
the actual amplifier. The requirements are that the range 2A be sufficiently 
small so that the characteristic is nearly linear within this range, and that the 

equilibrium voltage must be large compared to M so that the characteristic is 
almost linear up to the linear range of f. (This last requirement is not too 
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INTEGRATOR 7, “EXPONENTIAL STEP 


6. (t)=™M 


important; it will be obvious to the reader as he reads on how the nonlinearity 
could be taken care of.) 


It is desirable to put the problem in nondimensional form. Let 


t = AE 
Ar(E) = u(t) G3) 
yA, PA 
FE The 
Then 
1, r(E1 
+ flat)] = > ePrr(E)] = fr(e)] = { r(2), [r(@)|<1 
“1, r(E)<-1 
Hence, inserting Eqs. (2) and (3) into (1), one obtains 
rE + 1) Sar h{r(x)] dx, —& >0 a 


r(E) = 5 | Chat ent 


THE SOLUTION OF THE LOOP EQUATION 


From Fig. 3 it is obvious that the response of the system breaks up natu- 
rally into two parts. The first part is the slewing mode where the error, u(t), 
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decreases linearly from M to A, after an initial time lagA, or in the nondimen- 
sional formulation the error, r(&), decreases linearly from M/A to unity after an 
initial time lag of unity. The second part starts after the error signal has 
entered the linear region of f or h. 


The first part is easily calculated from Eq. (). One finds that 


¢ = MWA-1,y (5) 
1 p 


where Ey is the time at which the signal first enters the linear region, a 1 
corresponds to t) in Fig. 3. 


Because of the time delay, the error signal, r(§), continues to decrease 
linearly up to time fa £1 + 1. At Es €5: + 1 the error signal has penetrated 
the linear region of r by the amount. This is shown in Fig. h. 


Fig. 


We can now proceed to solve the second part. We initially restrict our- 
selves to the special case |r(§)| “1 for § = 6), i.e., once the error signal 
enters the linear region it does not leave it again. Redefining cur time origin 
so that § = 4) corresponds to T= O and r(6&) = w(t), Eq. (h) becomes 


wc +1) = (=p) - pf w(x) ax <>0 
7 o 


(6) 
wt) =1- yr O<tT<l 

One takes the Laplace transform of Eq. (6) using the boundary condition on w(t) 

for 0 < t<1. Taking the inversion integral one finds 


+ jo 


=] = 2: | 
w(t) Orr Jess 3(se° z hy 


g(r + iD: 
a ds, t >-0. (7) 


We now proceed to evaluate the above integral. One can show that the above 
integral is equivalent to a contour integral in the s-plane. The contour is a 
rectangle with corners at (0, -iNT™), (0, iNn), (-N, iNm), (-N, -iNT). As N goes 


\ 
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to infinity through odd positive integers the contour integral approaches the 
above integral (Appendix I). Hence, one can use the residue theorem. 


The zeros of the denominator are at s = 0 and the zeros of g(s) = se5 + pw. 


The pole at s = O simply cancels out the first term, 1, in Eq. (7). In order to 
find the zeros of g(s), let s = E + if; after some manipulation one finds 


~—& = n cot n, sin 4 # 0 (8) 
Eee = -Cos n, cos n # 0 (9) 


A few of these curves are plotted in Fig. 5 for various values of « and the 
zeros are circled. Of course, the families continue indefinitely in both the 
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positive imaginary and negative imaginary directions. The apparent zeros on the 
imaginary axis do not exist because of violations of the conditions on Eqs. (8) 

and (9). The zeros are all simple except for the special case #= el, & = -h, 

which is a double zero. One finds that the system is stable for 


n<n/2=1.57. 7 (10) 


For|n| + 2Tand for p sufficiently large, say 4 = 0.7, one can get very good 
approximations to the roots by linearizing Eqs. (8) and (9). 


Re(s.) = & = - In (2n+1/2) 7+ In yp (11) 


ny = (2 +2/2) m+ ee 


n=l], 2, 3 --- 


Im(s, ) 
(12) 


The principal zeros, i.e., those for which |q79| < 1/2, can be obtained by solving 
Eqs. (8) and (9) simultaneously by successive numerical approximations. A per- 
turbation expansion in powers of (1/2 - #!) is used to obtain the zero order 
approximation. Values for § o> No as functions of » are given in Table I and are 
plotted in Fig. 6. 
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Knowing the zeros of g(s) one can calculate the residues of the integral in 


rEqeel(7)< 
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» 1, 70 


(13) 


(Us) 


one finds three cases: 


Case ls: O<p<ert; two real zeros, call them aj and ao; all others complex: 


e e 
Se eee + 2 »y R 
w(t) = “SF a, Fl city (15) 


Case @: P= e-l; one double real root; all others complex: 


w(t) = 2e (r+ 1) (c +2) +2 z R. (16) 
n=l 


Case 3: e~l¢ y < 1/2; all roots complex: 
oo 
wit) = 2-2 ay (17) 
n=0 


The three cases correspond respectively to overdamping, critical damping 
and underdamping. We are interested primarily in Case 3. 


For the underdamped case, Gase 3, we can rewrite Eq. (17) in the form 


Etsy +1) cos Ca, aoc (4, - ,) J ee 
Jey eae cs) 


The error term, €(T), is bounded by (Appendix IT) 


w(t) = 2e 


+1 
lel <4 (2) aie oe (19) 


€(T) becomes small very rapidly; one sees that for the worst case, p= 1/2, that 
even at T= 2, |E(T)| < 0.06. 


The above analysis has been a linear analysis; it is valid only if |w(t)| <1 
for all1T>0O. In order to find the range of « for which the above is fulfilled 
one can go back to the defining integro-difference equation -= Eq. (6). This can 
be integrated piecewise, i.e., first for0 ¢ T<1, then for 1S T¥ 2, etc. 
Performing the indicated operations for 0<T< 2, then differentiating the 
resulting expressions for w(t), one finds that w(t) has a minimum at 


ae | f 2 
moe i 1<t,.S (20) 


and that, in this case, 
w(t) = 1/2 - p (21) 
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These equations show that |w(t,,)1< 1 for 4 < 1.5, and hence, that the linear 
analysis is valid for ak TA [Equations (18) and (19) show that w(t) does not 
go out of the linear range at any subsequent maximum or minimun.| 


We have thus shown that Eq. (18) is valid for all T> 0 for ¢1.5. Since, 
for stability, one must have h<7/2 = 1.57, the only range of interest for which 
the above analysis is invalid is 1.50< p< 1.57. This range is too smell. to be 
of engineering interest and we do not choose to consicer it. 


CONCLUSION 


Finally, we compute the minimum time necessary for the error to become 
essentially zero after starting out at an initial value M as a function of the 
parameter. Let "essentially zero" mean that w(t), as given by Eq. (18), must 
be less than or equal to Sfor allt >T,- We ignore the cosine contribution in 
computing T.. Then the time, T, in units of A, for the error to decrease to 
essentially zero from an initial error of MA, again in dimensionless units, is 
simply §], as given by Eq. (5) (i.e., the dimensionless time necessary for the 
error signal to get to the linear region) plus Te as given above. T is plotted 
as a function of » for 5 = 0.1 and 0.2 and for M/2A= 3, 6, 12, 25 in Fig. 6. 
The values are tabulated in Table I. The minimm can be obtained graphically or 
by examining the values of T. It is apparent from these curves that the minimums 
are shallow and therefore the value of p is not critical. 


The most important application of the results is not. the engineering one. 
It is well known that there are various approximation techniques which can be 
used to obtain engineering answers for the system analyzed here. They can also 
be used to engineer more complicated systems of roughly the same nature. But 
the importance of the solution obtained is that it can be used as a check point 
for current approximation techniques and others which will be developed. The 
utility is enhanced by the simplicity of the approximate solution and the fact 
that it was possible to obtain a simple quantitative bound for the error. 


APPENDIX I 
Consider the integral 


fle tate) | e(t +1) gg 
N 


s(se* + 1) 


where Ty is the rectangular path in the s-plane with vertices at (0,-iNT), 
(0,iNnr), (-N,iNtT) and (-N,-iNT) and which is suitably indented to enclose the 


origin. Let 
SNe | PeNiN | ENGL © fANr. 
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It is obvious that J] approaches the integral we want to evaluate as’ 
N+0o « 


Now consider Joe 


-N+iny | .s(t + 1)| 


J |ds 
| etree 3 


|si|se" + wf 
Let s = X + iNT, then 


[se + »| > e°|x sin Ny + Nr cos Mr | 


after dropping the imaginary part. Letting N be an-odd (positive) integer, this 
becomes 


Jee? + pl > er 
and thus 
x(t aids 
Nie ete 
Wal $ Jo (in) Nee aa 


Therefore, for T>0O 


N-+>co 
J), behaves in éxactly the same manner. 


Finally, consider J3, 


ite Biieleareo or: 
SN Rees . 


N+iNn | s| | se” + nl 


Let s = -N + iy, then 
|se~ + nl > lp - e Nin cos y + y sin y || 
after dropping the real part. Since |y| < Nr 
| eo Nin cosy +y siny| < on N(n + 1) 
Given a positive number Stuy there exists N, such that for N> Np, 


oNN(r +1) <n- 8 
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Hence, 
|se® +p] > 8 >0 


and therefore, 


Nr .-N(t + 1) 
a 2 -N(t + 1) 
[5,1 < L Saag ee dy = a e si oaNe 


Therefore, for T20 


N+>co 
APPENDIX ITI 


e(t) = 2 z eat oy eee meen Cac + (a - %) 


ire Sa + 1)* + nn 


Since one can satisfy himself that An oily Tl = als loieece, DY Looking atric. 5, 
then certainly 


(< +1) 
Sane 
n=l 


le(x)| <2 


From Eq. (11) 


co t+1 t++1 © c+1 
eos2 2 |e e(e) 2) 
n=1 (2n + 5)n 


This series can be shown® to be bounded by 


co tT +1 [s-e) 


1 dx a 
n=] ™ gt ti v 


Hence, 


le 
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THE STABILIZATION OF NONLINEAR SERVOMECHANISMS ENCOUNTERED 
IN ANTENNA INSTRUMENTATION* 


Jack Bacon 
Ohio State University 
Columbus, Ohio 


ae -- This paper investigates the problem of compensating the loop-gain of 
instrument servos to provide uniform transient response over a wide dynamic 
operating range. Particular consideration is given to applications in the realm 
of antenna measurements. A distinction is made between situations where the 
nonlinearity results from characteristics of the follow-up device and where it 
results from the loop-gain being functionally related to an external variable. 
In the former case the loop-gain is made self-linearizing; in the latter the gain 
is effectively equalized with an auxiliary logarithmic servo. This necessarily 
requires a suitable slave potentiometer in the nonlinear loop, actuated by the 
logarithmic servo. Certain advantages resulting from the use of ladder attenu- 
ator as the slave component are described. It is shown how this combination 
logarithmic servo and ladder slave attenuator can make a loop-gain invariable 
which otherwise would be functionally related to a variable E by the equation 

G = E*N, The idea is extended to the generation of polynominal terms of the same 
form as G. 


INTRODUCTION 


Automatic instruments employed in making antenna measurements are princi- 
pally those which record either logarithmic or square-root functions of amplitude 
or linear functions of phase. The major design difficulty encountered in these 
instruments is that of providing for uniform servo response throughout the full 
dynamic range. In the usual linear servo this consideration is no problem 
because the loop-gain is constant. In antenna instruments, however, the loop- 
gain is often a function of a generalized variable u. 


Figure 1 shows a typical situation where the follow-up device is function- 
ally related to the variable u. For certain types of amplitude plotters u is 
identified with the outputeshaft position. In other cases, say for phase 
plotters, u is identified with the input-signal level. This makes u independent 


FOLLOW - UP 
COMPONENT 
GAIN= G(u) 


Fig. 1 - Nonlinear SerVO0 


*The research reported in this paper was sponsored in part by the Air Research 
and Development Command, Wright Air Development Center. 
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of the output~shaft position, since the output is dependent upon phase of the 
input signal but not upon intensity. . Cases can be assumed where the variable u 
is completely divorced from the servo altogether. Since G(u) is a multiplier in 
the servo loop-gain, a change in this quantity can be expected to cause serious 
deterioration of the servo performance, independently of how u is related to the 
servo itself. Under actual operating conditions, variations between the extremes 
of complete insensitivity and a state of self-oscillation are usual. This paper 
is intended to show in a practical way how the influence of u on the servo loop 
can be nullified. 


SOURCES OF NONLINEARITY 
Consider an amplitude plotter where the output is proportional to a nonlin- 


ear function of input and the follow-up component is a potentiometer having fixed 
excitation. The pertinent example shown in Fig. 2 is a square-rooting device. 


Fig. 2 ~ Uncompensated square-root servo. 


Assuming a perfect servo, the functional relationship Eo = k@* yields Square-root 
response. Small-signal gain of the follow-up component, G(u), is defined as 
dE2/d6. In this particular instance G(u) = K@, which shows the parameter u is 
identified with 9. 


Servomechanisms using nonlinear potentiometers in the manner shown may have 
an incremental gain which is not simply proportional to the Slope of the follow- 
up potentiometer function. Certain other factors must often be taken into con- 
sideration in determining G(u = 6). For instance, the’ impedance of the input 
generator may be a contributing factor in the determination of G(u). It is 
worthy of note that under certain circumstances even a so-called linear servo 
may exhibit nonuniform incremental loop-gain. It is not within the framework of 
this paper-to consider how the various factors influence G(u) but rather how the 
servo may be made to have uniform dynamic stability in the presence of such an 
element. This discussion does not imply that the incremental loop-gain is always 
a variable in a servo designed to plot a nonlinear function of input-signal 
amplitude. Consider the particular type of logarithmic recorder shown in 
Fig. 3.122 Here the input voltage Ej is applied to an exponential attemator 
through amplifier K. The pick-off voltage e is servo-positioned for a constant 
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Fig. ) - Rudimentary phase servo. 


signal level so that the rectified voltage e' equals the reference voltage Vref. 
This makes Eqs. (1) and (2) hold. 


e = KE; € aoe =. constant (1) 
6 = C log (Kt E;) (2) 


Symbols C and K! are constants. The incremental gain of the attenuator is given, 
as before, as the derivative of Eq. (1) with respect to ©. This becomes a new | 
constant. Thus, an exponential attenuator used in this manner behaves for small 
changes in signal level as if it were linear. The servo performance is conse- 
quently independent of input-signal level over the usable dynamic range of the 
instrument. Recorders constructed using this principle confirm this fact. 


There are cases where the independent variable u is identified with the 
input-signal intensity but not with the sid Reece position. A servo designed 
to plot rf phase is a particular example.3»4 In an elementary design shown in 
Fig. h, a balanced hybrid junction is employed as the rf-phase discriminator. 
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Under proper operating conditions the magnitude of the audio error voltage & for 
small values of 6 is given by Eq. (3). 


3 = KE) Ed (3) 


The angle 6 is the amount by which the phase of Ey and Eo differs from quadra-~ 
ture. Thus, with a phase difference of 7/2 radians, & = O regardless of the 
amplitude of Ej. The instrument operates by using § to actuate a motor, which 
in turn moves an rf-phase shifter so as to minimize this error. Thus, as the 
phase of Ej varies, the rf-phase shifter "tracks" to maintain Eo at a constant 
phase difference of 7/2 radians. The phase-shifter position is. taken as an 
indication of the relative input phase, by actuating a chart mechanism from the 
mechanical output. Since the rf phase is a linear function of the shaft posi- 
tion, the static sensitivity of the hybrid junction is proportional to d_ &. 

dé 
The reference voltage Eo is normally constant in value, and, consequently, the 
gain function in this case is proportional to the magnitude of Ej, as is shown in 
Eq. (4). Thus, u is identified with the input-signal strength and independent of 
the shaft position. 


G(u) = 42 =¢ a = KE} (4) 


Inasmuch as the magnitude of Ej may vary widely in actual operation, the utility 
of the instrument in this simple form is necessarily restricted. 


COMPENSATION FOR NONUNIFORM GAIN 


Thus far we have discussed how servo loop-gain in some pertinent examples 
of antenna instrumentation can become a function of an external variable. At 
this point it would be well to consider some remedial measures which can be 
taken to maintain uniform transient response over the dynamic range of a servo 
characterized by including nonlinear elements of the type described.5 Figure 5 


AMPLIFIER 


COMPENSATOR 
GAIN = C(u)=G'(u) 


EXTERNAL CONTROLLING 
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NONLINEAR FOLLOW-UP 
COMPONENT GAIN Glu) 


Fig. 5 = Compensated servo with externally 
controlled nonlinearity. 


is a modification of the servo shown in Fig. 1. Under static conditions, where 
the input is maintained constant, Eq. (5) holds. It is clear in this case that 
the loop-gain is proportional to the quantity C times G. Thus, for small signals 
the effect of the variable u on the system 
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d pe aS ea : 
as E se G(u) (5) 


can be nullified if the quantities C and G are reciprocals. 


The first of two separate cases to be considered is where the independent 
variable u is identified with the output-shaft position. The second is where 
the independent variable u is an external voltage source. 


When the variable u is identified with the output-shaft position, the com- 
pensator C can be simply an additional potentiometer attached to the output 
shaft. Consider again the square-root servo of Fig. 2.° It was pointed out 
in this example that G(u = @) is proportional to 6. To obtain a reciprocal 
relationship between C and G, it mst be true that C = K/6. Figure 6 shows the 
compensator C attached to the output shaft. The reciprocal gain function is 
approximated by a set of shunts attached at judicious points along the 
resistance. 


It is obvious that a square-rooting servo can never be made to plot pre- 
cisely to zero, because the follow-up gain G(€) becomes vanishingly small in the 
neighborhood of 6 = 0. To maintain constant loop-gain in the neighborhood of the 
origin would require that the compensator gain increase without bound. Certain 
compromises can be made in this region, however, without seriously affecting the 
accuracy of the instrument. For example, modification of the law of response by 
altering G(@) so that it has a finite slope near the origin can be used without 
causing serious error. Since in this region there is usually noise present which 
tends to obscure the true reading, a slight distortion of the law of response 
here has no significant harmful effect. An aiternative in this situation is to 
stop the instrument short of zero. We have done this with instruments con- 
structed at the Ohio State University Antenna Laboratory. 


Uniformity of transient response across the dynamic range of the instrument 
depends upon the accuracy with which the compensator is adjusted. Figure 7 shows 
the uniformity of response which can be achieved over an 80 db range of input 
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signal or a lO db range of output. The response has proved to be sufficiently 
uniform so that there is little noticeable difference between its behavior anc 
that of a so-called linear servo. 


This same technique has also been applied to an instrument servo where the 
loop-gain. varied_exponentially with output-shaft position over an amplitude 
range of 100 db.‘ The compensator consisted of another exponential attemuator 
attached to the output shaft and having a gain function which varied in an 
inverse manner to the gain of the follow-up component. The finished design 
exhibited a uniformity of response very little different from that of a truly 
linear servo. 


In situations where the loop-gain is functionally related to an external 
variable, say a voltage, the problem of stabilization becomes somewhat more 
difficult. <A practical sclution, which places no severe constraints upon changes 
in the variable u, consists of employing an auxiliary servomechanism to obtain 
the necessary compensation. A suitable slave potentiometer mounted on the output 
shaft of the auxiliary servo can provide the required compensation to the princi- 
pal servo loop. Harmful effects, from the inclusion of a time dependent variable 
(viz., the compensating potentiometer) into the principal servo loop, can be 
satisfactorily minimized by making the compensating servo response relatively 
fast in comparison to that of the principal servo. This property can be readily 
designed into the system. 


The phase plotter shown in Fig. is an example of where an auxiliary servo 
can effectively provide the required loop-gain stabilization. The necessary 
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Fige 9 - Compensated phase 
Fig. 8 - Compensated phase servo. servo response. 
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modifications are shown in Fig. 8. It was previously pointed out that the error 
comparator has a gain proportional! to E; but independent of 6. As a result, the 
-compensator-gain function must necessarily be of the form K/E, throughout tre 
dynamic range of the instrument. The auxiliary servo is designed to have this 
property. Figure 9 shows the actual dynamic response of an instrument employing 
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compensation of the type discussed and having an rf aseke amplitude free to vary 
over a lO db range while the phase is being recorded. 


SERVO COMPENSATOR DESIGN 


Technical investigations often lead to original by-product results with 
usefulness extending beyond the initial scope of the inquiry. Such is the case 
here. The results to be described grew out of a need for a suitable compensating 
servomechanism for use when the independent variable u was an external voltage. 
In viewing the desirable properties of a suitable compensator, the following 
points are salient: 


1. It should employ servo techniques. 
2. It should be capable of operating with an audio input voltage. 


siege evace be capable of operating over a wide dynamic range (80 db of 
audio). 


. It should be as flexible as.possible so that the same basic design may 
be adapted to a wide variety of compensation requirements. 


5. The functional accuracy may be relaxed in comparison to that normally 
required of an accurate amplitude plotter. 


The requirement for a wide dynamic operating range immediately leads us to 
consider a modified logarithmic servomechanism for use as a compensator. The 
question which logically arises at this point is: what sort of a slave component 
will be suitable, and can one form of it fulfill a variety of requirements? 


Figure 10 shows a device which has great flexibility in the type of gain 
function it can simulate and, in addition, has other important features which 
will be emphasized. The top portion is simply an audio logarithmic servo wherein 
the follow-up potentiometer is exponential and has an attenuation ratea. The 
output shaft, having a normalized angle of rotation © (0< © <1), is coupled to 


MASTER ' POTENTIOMETER 


REFERENCE 
VOLTAGE 
B 


Po 
V 
eee 


Fig. 10 - Logarithmic servo compensator. 
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a sleve potentiometer. The slave is also an exponentially tapered attenuator 
but has in general a different rate of attenuation than the master potentiometer. 
Let us define its attenuation rate as Na. The variable u is represented by the 
voltage E, applied to the master potentiometer. The servo positions the pickoff 
on the potentiometer so as to make its voltage ¢ constant. Thus, 


ee By €%°© « constant - (6) 


The input signal to the compensator is the voltage V which may be applied to 
either end of the slave potentiometer. Assume that it is first applied to the 

A end. The output voltage from the compensator is the voltage v taken from the 
tap O on the slave potentiometer. Since the voltage distribution along this 
potentiometer is also exponential and has an attenuation rate Na, the compensator 
output v is related to the input V by Eq. (7). 


vr Viegas (7) 


The compensator-gain function C(E]) is the ratio v/V. Combining Eqs. (6) and 
(7) and rearranging term gives for C(E}) 


C(E,) = v/V = Key ™ (8) 


where K is a constant. Thus, with this technique we are able to generate a gain 
function proportional to a voltage E; to the minus N power, the N being deter- 
mined by the ratio of the attenuation rates in the two exponential attenuators. 
No constraints are placed on N except that it be a positive number. 


Now consider a slight modification.to the circuit. Transfer the compen- 
sator-input V from terminal A on the slave attenuator to terminal B. Equation 
(9) is valid in this case. 


Wi Vig -Na(1-6) (9) 


Again substituting from Eq. (6) and reordering the terms yields 


C(Ey) =v = KE," ° (10) 


Thus, the method described will lend itself to the generation of gain functions 
of the form given in Eq. (11). 


C(u) = mu (11) 


Here u is the generalized control voltage previously discussed in regard to the 
compensation of loop-gain. It is clear that a wide variety of compensation 
requirements can be fulfilled using this technique and that more complicated 
compensation functions could be synthesized by combining the output of several 
slave components where the exponents in general were different. 


The reader will recall that in the phase plotter it was required that the 
compensator gain be proportional to K/u. This was obtained by constructing both 
attenuators to have the same attenuation and by using terminal A and the arm 
respectively as the input and output of the compensator. The logarithmic prop- 
erty of the master servo made it possible to obtain easily an 80 db dynamic range 
on the input control voltage without any difficulty. 
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An obvious extension to the described technique is that of using the compen- 
Sator as a function generator in its own right. This can be done by applying a 
fixed potential to either terminal A or B of the slave unit and using the pickoff 
voltage as the output, thereby synthesizing functions of the form given in Eq. (12) 

tn 
v= KEy (12) 
It should be possible, in principle at least, to construct polynominals of the 
form shown in Eq. (13) if a number of the slave units are employed, each 


joe em 
v=). Apo (13) 
j = =n 
having a different rate of attenuation. This would, of course, require inde- 


pendent amplifiers to adjust the separate constants and an adder to perform the 
summation. 


A further point worthy of note is the flexibility of the design when used 
either as a function generator or aS a compensator. Recall that in the slave 
potentiometer the exponent of € is N46. Here N can be associated with either 
Nor 6 By associating N with 6, the law of response can be altered merely by 
using a gear ratio between the master potentiometer and the slave instead of 
direct coupling as shown in the figure. 


CONCLUSION 


In review, we have discussed how certain nonlinearities naturally arise in 
some important antenna=-instrumentation problems. The nonlinearities were divided 
into two categories. First, there weregthose which were functions of the output- 
shaft position; then, those which were dependent upon an external potential, 
which could be independent of the controlling voltage of the servo. In the first 
case, it was shown that an additional compensating potentiometer could be affixed 
to the output shaft and coupled into the error path of the servo to bring about 
satisfactory gain stabilization, even though the principle was based on small- 
signal response. In the second case, we saw how an auxiliary compensating servo 
could be employed successfully to compensate for external influence on the non- 
linear servo. This required a technique for obtaining a wide variety of gain 
functions to meet diversified conditions. The compensator was designed to obtain 
the desired compensation under a variety of conditions. It was indicated how a 
slight modification of the compensator would convert it into a function generator 
for the purpose of generating terms of a constant-coefficient polynominal. There 
is much work remaining to be done before the full possibilities of the function- 
generator technique can be effectively assessed. However, at the present time it 
offers some challenging possibilities. 
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ON THE DESIGN OF A-C NETWORKS 
FOR SERVO COMPENSATION 


Harold Levenstein 
The W. Le Maxson Corporation 
New York, N. Y. 


Summary -- This paper presents an analytical method for analysis and synthesis of 
networks for AC servo compensation, 


The response of a linear network to a modulated suppressed=carrier excita- 
tion is formulated in terms of in-phase and quadrature carrier components. The 
relationship between the modulation on these components end the exciting modula- 
tion is shown to depend upon operators simply related to the original network 
function. 


The expansion of the data-frequency operators into partial fractions is 
shown to lead to a simple synthesis procedure for deriving the original network 
operator from the in-phase or quadrature operators. 


As an example, the process is applied to the derivation of a represent= 
ative network for lead compensation of an AC servo responsive to the in-phase 
component of error signal. 


The use of RC networks in AC servo compensation is shown to be limited to, 
derivative types of equalization. 


INTRODUCTION 


The carrier-frequency servomechanism occupies an increasingly important 
place in the instrumentation field. Evidence of this appears continually in 
advertisements and in papers describing developments. A continuing problem 
in this field is in improving methods for analyzing end selecting networks in 
order to stabilize or compensate carrier servos without resorting to demodula=- 
tion techniques. Starting with Sobczyk's work in 1945, the emphasis has been 
on frequency analysis techniques, end methods of approximation based upon fre- 
quency response.» 2 Probably the most popular end effective method tc date 
has been the use of low-pass to band-pass network transformations as a means of 
finding appropriate networks. 


In this paper we teke a different approach. Using Laplace transform 
techniques we first develop exact relations for network response to suppressed 
carrier modulation, The resulting relationship is examined in terms of the 
pole-residue paremeters of the network. From this we develop specificetions 
that a required modulation operator must meet in order to be realizable. 

Using the pole-residue relationships, we develop a method for modifying a given 
operator in order to make it realizable; a technique for recovering the realiz- 
able data frequency operator is then described. 


The results are applied to the analysis of the general performance of RC 
networks and to derivation of a network transfer function providing lead-lag 
operation on the modulation. 
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The work assumes that a servo motor is equivalent to a demodulator re=- 
sponsive to the modulation on the in-phase component of the carrier. 


THE DATA FREQUENCY TRANSFER FUNCTIONS 


It is assumed thet the input signal is a suppressed cerrier modulated 
signel, as commonly encountered in AC servos. It may arise from the error 
signal of a synchro control transformer or as the consequence of mixing com= 
mand voltage and feedback voltage. The essential feature is that it may be 
represented mathematically as 


@g = e(t) cos wot (x) 


The situation is depicted in Figure l. 


Fig. 1 - Diagram of problem. 


In Appendix A we derive the following result for the output of the 
networks 


@o= Cos Wot ES: E(u) exo(ut)4 2 (wu + iw,)+N(u - =))\ 
(2) 


+ sin vot : [Ls du E(u) exp(ut) {3 (N(u + iw.) — N(u - ama 


TT i: 


In this equetion E(u) is the Laplace transform of e(t), and N(u+iwg) 
and N(u - iw,) are the transfer functions of the network evaluated at 


p=(ut+iw,) and p=(u - iw,). It appears from this expression that the out- 
put of the network can be regarded as the sum of two independently modulated 
suppressed=carrier components, one in-phase with the original carrier, and 
the other, the sin term, in quadrature with the carrier. The modulation on 
each component is seen to be the inverse Laplace transform of E(u) multiplied 
by a suitable function of u that depends on the network. We may regard these 
functions as the date-frequency operators of the network, An equivalent cir- 
cuit is AN in Figure 2. 


This formulation is quite general. It has not been necessary to specify 
that the response is a current or a voltage, so that N may represent a driving 
point impedance or admittance, or a transfer function. It merely states that 
if the system function relating two variables in a linear Re tayoris is known, and 
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e(t) coswit 


Fig. 2 - Equivalent data frequency operation. 


if one of the variables is suppressed-carrier modulated, the other consists of 
two suppressed=-carrier modulated components, whose modulation is uniquely re- 
lated to the reference modulation through the pair of operators. For conven- 
ience we define the operators as 


d= [ N(u+iwe)+N (uv - iw,)| | (3) 


Dg=S [ w(u+ iwc) -N (u = iwe)| (4) 


where D; and Dy are the in-phase and quadrature operators respectively. 
ADDITION AND MULTIPLICATION PROPERTIES OF THE OPERATORS 


If two networks are connected so that their responses to an input 
signal are added to form the total response, the data~frequency operators 
add to form the net data-frequency operators. 


ea ero De 


i: (5) 


However, if two networks are connected in cascade, their data-fre- 
quency operators do not multiply, as do the original network system functions. — 
The overall system function must be calculated in order to determine the over- 


all data-frequency functions. 


Di= $ [ma (ating) No(utiw,) +N, (u - iwg) No(u - ing)| 
: | (6), 
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Fortunately, the addition and multiplication operations are commutative. 
That is, the order in which the networks are added, or the order in which they 
are cascaded, do not alter the overall data-frequency operators. 


It will appear later that these properties dictate the manner in which 
the network properties are synthesized. 


PROPERTIES OF THE DATA-FREQUENCY TRANSFER FUNCTION 


The data-frequency operators are the sum and difference of N(u+iw,) and 
N(u = iwc). . 


These functions are complex rational functions of the complex variable 
ue The original transfer function N is a real rational function of p. It may 
be written as 


~~ 


P(p 


7 
Q(p i 


| 


N(p) = 


~~ 


where P and Q are real polynomials in p. 


The roots of Q(p) (the poles of N) must lie in the left half of the 
p plane. Real roots lie on the negative p axis; complex roots occur in 
conjugate pairs. 


The roots of P may lie anywhere in the p plane, with the proviso that 
complex roots occur in conjugate pairs. 


The degree of the polynomial P must be equal to or less than the de- 
gree of the polynomial Q. There is also a constant multiplier associated 
with Ne. The maximum permissible value depends upon the network topology by 
which N is to be realized. For the purposes of this analysis we may ignore 
the multiplier. However, in the determination of network component values, 
it would be necessary to consider the multiplier. 


Under the transformations p=u+iw, and p=u - iw,, the poles and 
zeros of N are shifted down and up respectively in the complex u plane. 
This situation is shown in Figure 3, It is obvious that for every root of 
N(u+iw,) there is a corresponding conjugate root in N(u = iw). 


If the zeros of N are denoted by Greek subscripts and the poles by 
Latin subscripts, we have the following relations: 


N(p) = (2 - Pap - pg) sors 


(Pp - Pa (Pp = Pp ) eeee (8) 
S i (u ne Pot ive) (u ae Pg tiWc) ore 
ae 2|{(u - Patiw,)(u = Pptiwy) ces 
(U = Dy- iwe)(u = Pg 7 iw.) esse 
(u = pam iw,)(u - Pp - iw,) ea (9) 


N 
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i|(u = pytiwe)(u - pgtiwe) ..6- 
= 
q 2 (u = Pat iw,)(u = Ppt+iw,) ecee 


(u = Pe tiw,)(u — Sepa iw.) e008 


(u ed Pa~ iw,)(u St Pp a iw,) ©ee.n 


(10) 


There are three distinguishable root locations. If a root of N is 
on the real axis, p=+a, it appears as a complex root u=a - iw, in 


N(u+iw,) and the conjugate u = a+iw, in N(u = iw,). 


If a general pair of complex conjugate roots p = (atib) appears in N, 
it transforms into a+i(b - =i and a - i(b+w,) in N(u - iwg), and a+i(b+wg) 
and a - i(b = wo) in N(u+iw,). ; 


pole on complex real pole 
+iWe line. pole 


Fig. 3 - Movement of typical poles of N. 
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Finally, if the complex pair takes the form p=(atiw,), they transform 
into a, a-2 iwg and a, a+2 iw,. That is, both N functions contain the factor 


(u- a)e 


Attention will now be concentrated on the ae function. 
Di=t | nlu + iw,)+N(u - iw.) 
i=%> N(u + iw,)+ N(u iW, 


aah P(u+ iwe)  P(u = iw.) 
2 1Q(u+t+ iwc) Q(u = iw,) 


fe P(u+ iw) Q(u - iwe)+P(u = ine) Q(ut ing) ; 
=— 11 
ee Q(u+ iw,) Q(u = iw) 

Reh 

— Qi 


Q; contains all the roots of Q repeated twice, once shifted up iw, and 
once shifted down iwge Complex roots occur in conjugate pairs. Therefore, _ 
Q3; is a real function of ue. Since the real parts of its roots have not been 
altered, all of the roots lie in the left half u plane. If a root transforms 
into a real root it appears once in P; and twice in Q;. Thus, the degree of 
the polynomial Q; is equal to or less than twice the degree of Q. It is exactly 
equal to twice the degree of Q less the number of roots of Q having an imaginary 
part equal to w,. A real root in Q transforms into a complex conjugate pair. A 
complex conjugate pair transforms into a pair of complex conjugate pairs, except 
for the special case when it transforms into a triple consisting of a real root 
and a complex conjuzete pair. 


The roots of Pj are not readily established. Pj is the sum of two poly- 
nomials whose coefficients are complex conjugates of each other, so that Pj is 
a real function of u. 


If a root of P transforms into a real root, this becomes a root of Py and 
the corresponding term may be factored out. However, the roots may lie anywhere 
in the complex u plane without any simple relationship to the roots of Pj. The 
degree of the polynomial Pj is equal to or less than the sum of the degrees of P 
and Q on reduction t6 simplest form. 


PARTIAL FRACTION EXPANSION 


In determining the desired data-frequency trensfer function for a servo, 
it is customary to express it as a ratio of numerator factors divided by de-~ 
nominator factors. That is, ordinarily the poles and zeros of the desired 
operator are specified. While this determines a unique function, it is not - 
the only way to express the operator mathematically. This is fortunate, for 
as was noted above the zeros of the data-frequency operators are not readily 
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ascertained. For this purpose the partial fraction expansion is most useful. 


Instead of expressing N as a ratio of two polynomials, it proves desirable to 
express it in partial fractions, 


Ne= = 24 a + a2 a3 


=Pl p> pt ao (13) 


In this form the P; are the poles of N, a, is the residue at infinity 
and may be zero (when the degree of P is less than the degree of Q), and the 
a, are the residues at the poles of N. The a, generally are complex, occur- 
ing in conjugate pairs associated with the complex conjugate poles. The real 
parts of the a, may be negative. . 


Forming the expression for Dj in these terms 


a a 
N(u+iw,) =a + atalag + Pine tla ant + (14) 
ue Pi+ 1Wo ue- Po+1W, ece 
ay i) 
N(u - iw,) = %o+ ear aren (15) 
ee ogee ed Poe 
B | al al 1 
,= 4 fe me ee TI A ES ETRE comm acta 
Bir oe 2 Sercepeerr rn eenerertee a ah (16) 


If p; is real, then the typical term d 


eee ( a = sal | ) 
BIg « P}+iW, ei Sh ies iW, 


consists of a pair of complex conjugate poles. 


If py is complex its conjugate will also be found and the group of terms 
associated with it becomes 


a) ay BY att 
x = =, * 
u-= Py + iW, ss ial Sg iw, us Py S4iWoe UUs + iw, 


3 
qd =— (17) 
2 


This shows that in the Dy function every complex pole is accompanied 


by @ complex conjugate pole and the residues at these poles are also com- 
plex conjugate, derivable directly from the corresponding residues of N. 
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Thus, the essential transformation properties of Dj are these: the 
poles of N are carried over in duplicate (as described earlier) and the 
residues at the poles are carried over intact. The problem of inverting 
the transformation can be accomplished by means of them. 


THE INVERSION PROCESS 


The process of inversion, of finding N from D;, is primarily one of 
identifying poles and residues with the poles and residues of N(u+iw,) and 
N(u - iwg) and then shifting the poles as required. The rules are relatively 
simple. 


1. Express Dy in partial fractions. In the equation below, real and 
imaginary parts of the poles and residues are explicitly denoted in order 
to indicate the difference between real and complex poles. If Dj is a proper 
function, it will take the following most general form: 


y ay 


Da a5 a= ai F ing wat, > im, ces (18) 


terms corresponding to real roots in N. 


eotene 282 ao - ibs 
+ ——__—__ + ——__ + —_2_ 2 
u-Q+2 iw, ou = O u-@ - 2 iw 


terms corresponding to complex pairs (Lot iwg) 
a hay eaten Silas Mert as fer, EPs 


terms corresponding to general complex pairs 


Ast By (3s Peas 
) Veith 6 eee 


( o,, ti 
| 2 2 ; : 


2. From this group identify the partial fraction expansion of 
N(u+iw,) as follows: 


287 


arin (19) 


N(u+ iw.) =aot+ 


2(ag+iby) 2(a, = ib) 
+ + 


u - 5+ 2iw, u ~X, 
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2(ag+ ibs) 2(az = ibs) 
+ ———————— +, 
u “Az ti fez U - 0, “ ifg, 


3. Make the substitution u=(p = iw,). Then N(p) is given by 


we 2a} 
De hea ot rOeeeiet* (20) 


2(ao+ ib) 2(a5 = ibs) 
et 


ood @ 


p-AHp+iw, p-o = iw, 


9 2(ag +ibs) r 2(a3 = ibs) 
Pp -X.+ i(B, oad Wo) Pp -&, as i( 2,4 We) 


THE QUADRATURE OPERATOR 


The general properties of the quadrature operator are the same as 
those of the in-phase operator, in terms of poles and zeros. However, 
due to the i multiplier and the negative sign associated with N(u - iw,)>» 
the emphasis in the partial fraction expansion is on the imaginary part 
of the residues. The typical form for Da corresponding to equation 18 is 
displayed in equation 21. 


2 SS ee een 21 
Dq u- 1+ iw, u-d, - iw, (21) 


u “A+ 2iw, u -A> U =O, - 2iw, 
e bz - iag a bg+iag " bg+ ia, bz - jag 


weXsti fA, u-Os- ifs u-Gstia, u-t,- 17, 


Note that the constant term has disappeared and that the residue 
associated with the real poles is the imaginary part of the residue in N. 
The inversion process can be performed through equation 21. The constant 


part of N is then arbitrary. 
| AT 


AN EXAMPLE 


An example may serve to clarify some of the problems and illustrate 
the technique of inversion. 


Consider the following problem. The power phase of a two-phase 
servo motor is connected so that the motor is responsive to the in-phase 
carrier component, This is a common situation. It has been determined 
that the transfer function’ operating on the data should approximate that 
of the ordinary lead-lag network. That is 


ut a 
Dj; 2&—— where a<b (22) 


ut+b 


We append to this the requirement that Dg should be small at modu=- 
lation frequencies. This is desirable since the quadrature component acts 
to reduce the dynamic range of the amplifiers and to heat the motor with= 
out producing torque. 


Let 
uta 
Dis 3 adb (23) 
utb 
~ bea 
Di= oe! u+ b (24) 


It is clear from the previous discussion that this does not have the 
proper form, since every real pole must be accompanied by a pair of complex 
conjugate poles. Therefore, we modify Dj by adding to it additional terms. 


span bda rah | JUSS) FARE a | Stee eet 


utb 2 |u+b+2iw,| 2 |utb - 2 iw, (25) 


Ais arbitrary. We may use it as a means of adjusting Dg 

We can make a further adjustment to compensate for the introduction 
of the two complex poles. In the vicinity of u=0 (the spectral range of 
servo modulation) the two complex poles behave like constants. They are 
ordinarily located so far away from zero that we can add to the co.istant 
term an additional term designed to eliminate this zero frequency behavior 
and expect this situation to hold over a wide range. In this case, the 
constant would be 


2b (b - a)+4w.f 


1/2 


re A 
oe ae 
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Since our technique is not affected by this, the modification will not be 
introduced here, 


Od 
Dy is now in a suitable form for the inversion. 


N(u+iw,)=1 - ell eres el (b =a) = if 


(26) 
ut b+2iw, ued 
(b= a)+if (b- a) -14 
N(p) =1- Me SE Ser phaie: (27) 
pt We p+ iw, 
Calculate Do 
i : ; 
Dg= 5 | (uti) - N(u = iwe)| (28) 
_ i | -4iwe (b-a)+ i268 (u-b) 2 if | 
SEE RTL RES SENS eee an 29 
2 (u+b)?2+4 we utb ey 
2 Wo (bea) - A (utd) 2p 
is ST eee (30) 
(u+b)°+4 we u+ b 
The selection of es zero reduces Dg to 
2 Wo (b = a) 
LE aT RE TT (31) 
(u+b)* + 4 Ww 
a) bea 
= for uz0. (32) 
2 Wo 


For this selection of 2 » the modulation on the quadrature component 


is proportional to the input modulation, to (b = a) end inversely proportion- - 
al to the carrier frequency. 


From (27) with @=0 


(b- 2) (b= 4) 
N(p)=1- Spying ped — iW, (33) 
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p+ 2 apt (we+2 ab = vb”) 


(34) 
(p+ b)?+ we 


N(p) = 


USE OF RC NETWORKS 


It’ is sometimes a self-imposed requirement of the designer that he 
limit his selection of components for networks to resistors end capacitors. 
The data-frequency operators of RC networks necessarily have their poles 
distributed along the tjwc lines in the u plane where, as poles, they con- 
tribute little to the operations on the modulation spectrum. However, under 
these circumstances it is necessary to scrutinize the location of the zeros 
of the data-frequency operators, since by suitable selection of residues it 
is possible to make them the effective part of the overall operator. A case 
in point.is the parallel Tee network often used to provide data-frequency 
lead. The transfer function of the network can be taken as: 


p+ we 
p’+2 aWo+pt we 


N= (35) 


The numerator has purely imaginary roots at aa iw,e As required in an 
RC circuit the denominator has real roots, equal when a is unity. The fre- 
quency-shifting transformations carry N into the two conjugate functions 
N(ut+iw,) and N(u = iw,). 


u (ut2 iw,) 


Al a ee tone oF Bed Dee IS a Dae 
Mat awe) = Lewe (atVe® 14a) [ utwy (a- Ve? - 1+ 4)| 


u (u = 2 iw) 


[ w+ we (aty a2 -l- i [ utwe (a-V a -l- i)| 


ue +2 aWo u-+ 4 wut 4a we 
D=u (36) 
: ut+ Aaw, w+ u-(4 a@ we + we) +8 awe ut+4 a& we 
wi uif w.>u 3 
BoM cm (37) 


A better approximation in the low-frequency range is 
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1 
i Re can) Us Wo> u (38) 


This indicates the presence of low-frequency lag. Note that the scale 
factor associated with the differentiation operation is inversely proportional 
to the carrier frequency. This is not the case with the network designed pre- 
viously by way of example. There (equation 23), with a taken as zero, the 
scale factor is approximately 1/o. Since there is no point in controlling this 
break frequency as far out as b=wWe, it is possible to obtain higher gains with 
this circuit. Note also that the selection of B relatively small compared with 
Wo keeps the quadrature noise down. 


The quadrature function for the parallel Tee is 


~4aw vw 


he (39) 
q ut 4 4aw, u? + ue (4 a® we + we)+ 8 aw. u + 4 2 we 


Thus, the low-frequency behavior of Dq is approximated bys 


oe 4 2 
D.= -_ ra we u (40) 


LOW-FREQUENCY LAG CONTRIBUTION OF RC NETS 


In equation (38) it was observed that the denominator of the RC net- 
work data~frequency operator could be fairly well approximated by a simple 
first-degree function of u, contributing a lag term at low frequencies. 
Since all of the roots of the denominator are on the +jwe or =- jw, lines 
this. is a fairly good approximation. The time constant of this leg is the 
sum of the time constents of the denominator first-degree factors, that is, 
the sum of the reciprocals of the roots of the denominstor. 


We have then, T, the approximate time constant is 


T= -3 > (41) 


Since the roots x, are the real poles of N shifted +iw, we have 


T= nine 7 (42) 


where poles of N are given by (- Qj). 
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Thus, in the parallel T we have poles at p=W, (-2ty a’ - Te 


Substituted in (42) these lead to a time constant 


2 
aw re} 


T= 


Generally, if the designer does not require control over too wide a 
bandwidth, he can employ RC networks to generate any degree of derivative 
type compensation, eccompenied by a small, simple leg term. 


EFFECTS OF SHIFTING MOTOR REFERENCE PHASE 


If the motor reference phase is shifted through an angle 6 so that it 
is no longer wholly in quadrature with the signal reference carrier, the 
useful data-frequency operator is a linear combination of Dy and Dg. If the 
motor is responsive to the cos (wt +) component of the carrier the effective 
data operator is 


D = cos? Ds ~ sin O Das 


The effective quadrature operator is 


D= sin 9D; + cos 9 Dg. 

The properties of these functions are generally the same as those of 
Dye Since maintaining a specified phase shift different from zero or ninety 
degrees is usually difficult, it appears unlikely that much use can be made 
of this relationship. 


CONCLUSION 


It has been shown that the data-frequency operators of a network 
excited by a suppressed=carrier modulated signal are readily derivable from 
the network parameters. When a specification exists for such an operator, 
it is readily checked for realizability end various artifices may be em= 
ployed to modify it into appropriate form. Once in suitable form the pro- 
cess of converting the data-frequency operator into the basic carrier-fre~ 
quency operator is accomplished by simple algebraic means. There are usually 
at hend a sufficient number of arbitrary parameters to permit simultaneous 
minimization of the quadrature signal component. 


No special reference has been made to detailed methods of realizing 
‘the network transfer function as a network, since the literature has made 
many methods availeble, and the designer usually has his own preferred 
technique. (Reference 3) 


The use of RC networks: is shown to be directed primarily toward deriv= 
ative types of compensation. If control of lag operators is also desired, 
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it is necessary to resort to circuits possessing lightly damped poles, i.e., 
RLC circuits of high Q components. 


By means of the techniques described, it is expected that the variety 
of networks useful for AC servo compensation can be extended beyond those 
that have been exploited on empirical grounds. 


APPENDIX -- DERIVATION OF DATA=FREQUENCY OPERATORS 
Given a signal eg of the form 
es = e(t) cos wot. 


Assume it applied to a network whose system function is H(p). 
That is 


R(p) =H(p)' E(p) 


Where R(p) is the response, E(p) is the excitation, expressed as Laplace 
transforms. Then 


R(p)=H(p) L(eg) 
L(eg) = fot) cos W,(t) exp (=pt)dt 


rele 


fete) [ exp(-p -i wo) | dt 
+h foley [oxp(-p + 4 we)t] at 
=$ Ee: iw,)+E(p - ix,)| 


R(p) =5 Ee + iw.) H(p)+E(p - iw,) ap) 


Take the inverse Laplace transform 


j co 
r(t) = - Ee ihe E(p+iwg) H(p) exp(pt) dp 


ee 
tsar yf: ao? - iw.) H(p) exp(pt) an| 


In the first term above, let p=u - iw,; in the second let 
p= V+ iwc. 
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joo 
Jf, 30) H(u - iw.) exp (u - iw,)t du 


r(t)= = eh 
4.00 
+f xlv) H(v+ iwg) exp (v+ iwg)t ee] | 


Since the integrals are each definite, over the infinite range we may 
set v equal to u without altering the values of the integrals. Combining 


400 
r(t) =4f oe [30 H(u = iwg) exp(u = iw,)t 


+H(u + iw,) exp(u+ iw,)t 


du 


1 


1.00 
- i du E(u) exp(ut) Ee - iw,) exp( - iw,t) 
1 J/J-] co 
+H(u+ iw.) exp(+ in] | 
Let exp (iwct)=cos wot+i sin wot 


exp (-iw,t)=cos w,t - i sin wot 


Combine eid factor out the cosine and sine terms. 


r(t) = cos wot ee, ‘ an E(u) exp(ut). ; Ec iw,)+ B(uciwe)| 


+ sin w,t SREY a ace E(u) exp(ut). = Eee ‘iwg) -i(u-im)| 


For the transfer function of a network we have used the symbol N as a special 
case of He 
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FUNDAMENTAL EQUATIONS FOR THE APPLICATION OF | 
STATISTICAL TECHNIQUES TO FEEDBACK-CONTROL SYSTEMS* 


G. Ae Biernson 
Sylvania Electric Products Inc. 
Waltham, Massachusetts 


Summary -~ The basic equations necessary for applying statistical techniques to 
the design of feedback-control systems are presented. The autocorrelation func- 
tion of the output is computed by a transient technique which treats the input 
autocorrelation function as a transient input to the system.. By transforming 
this procedure equations are developed for relating the spectral densities of 
the input and output, and a means of performing the computation on an analog 
computer is presented. 


INTRODUCTION 


Scope 


The basic mathematical equations necessary for the application of statis- 
tical techniques to the design of feedback control systems are presented in a 
concise, yet reasonably simple manner. This paper gives a simple physical 
interpretation of the statistical calculations by showing that the various 
correlation functions can be related to each other by transient techniques. 


The analysis to be described contains a number of important theoretical 
gaps and hence is not mathematically rigorouse However, the purpose of the 
analysis is merely to develop the statistical tools necessary for engineering 
computation and to present a plausibility argument to give the engineer a physi- 
cal understanding of the concepts involved. The basic material of this paper is 
in no sense original. The author is presenting well-known statistical relations 
in a manner which, he believes, makes them more understandable and useful, to the 
engineer. For more detailed information the reader is referred to Truxal" and 
James, Nichols and Phillips.° 


This paper shows how to calculate the mean-square value of the response of a 
linear system from the autocorrelation function of the input and the system 
transfer function (or impulse response). It presents two methods for performing 
this calculation, a transient and a transform method. The transient method 
treats the autocorrelation function of the input as a transient input to the 
system and calculates the mean-square value of the output by determining the sys- 
tem response to that transient input. The transform method is derived by taking 
the transforms of the operations involved in the transient method. 


*The research reported in this document was made possible through the support 
extended the Massachusetts Institute of Technology, Servomechanisms Laboratory, 
by the United States Air Force (Weapons Guidance Laboratory, Wright Air Develop-— 
ment Center) under Contract No. AF33(616)-2038, Expenditure Order No. R5S5\-31]- 
SR-12, M.I.T. Project No. D.I.C. 7138. It is published for technical informa- 
tion only and does not represent recommendations or conclusions of the spon- 
soring agency. ' 
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Difference Between the Analysis of Noise and Signal 


Statistical techniques are very useful for studying the effect of noise in 
the reference input to a feedback control system. A detailed knowledge of the 
noise response is generally not necessary; usually a statistical measure is 
adequate. In fact, when a designer actually has a complete plot of the response 
to a noise input, he often condenses the data statistically in order to define it 
quantitatively. Consequently, it is often desirable to perform in a statistical 
manner the complete analysis of the effect of the noise by condensing the noise 
input into a convenient statistical form and from this compute the necessary 
statistical properties of the output. The statistical condensation of the input 
that is used is the input autocorrelation function 04;(T). Since this is much 
simpler in form than the actual input x;(t), it has much less information but 
does contain sufficient information for determining the mean-square value of the 
response of any linear system to that input. 


Although statistical analysis is usually adequate for determining the effect 
of a noise input to a feedback control system, it is generally not adequate for 
determining the effect of the signal portion of the reference input, which the 
system is designed to follow. One is usually interested in much more information 
about the system response to the signal portion than merely the mean-square 
error. Besides, the signal portion often cannot be considered statistically 
stationary, and a given piece of signal input data is often of such low frequency 
with respect to the period over which it is measured that a statistically mean- 
ingful autocorrelation function of the signal portion cannot be obtained. If an 
autocorrelation function is computed from the data available, it often does not 
yield reasonably accurate information concerning the mean-square error of the 
system response. On the other hand, a plot of the signal portion of the refer- 
ence input is generally fairly simple in form because the signal portion must be 
of low frequency with respect to the bandwidth of the system in order for the 
system to follow it. Consequently, there is no real need to condense the input 
information statistically before calculating its response. 


Thus, for feedback-control applications the statistical techniques should 
be used in general only for studying the response to noise inputs. When studying 
the signal portion of the reference input, the complete response can be calcu- 
lated quite readily by the techniques mentioned in an earlier paper which 
presents a graphical method for calculating the response of a system to an 
arbitrary input by expressing the response as a sum of error coefficient terms 
and transient terms. In fact, this method can also be used in calculating the 
statistical properties of the response to the noise portion of the input because, 
as will be shown, the statistical calculations may be considered to represent the 
solution of a transient problem. 


Statistical techniques have received a wider application in the field of 
communications than in the field of feedback control. One reason for this is 
that in communications the problems encountered have more of a statistical 
nature. It is often desirable to treat in a statistical manner the signal por- 
tion of the input to a communication system as well as the noise portion, because 
the signal portion generally has such a complex nature that it is very difficult 
to work with a complete plot and it is usually of such high frequency with 
respect to the period of measurement that an autocorrelation of a sample can be 


‘quite meaningful. Besides, communication applications generally require that the 


57 


output follow the input with far less accuracy than do feedback-control applica- 
tions, and consequently it is often sufficient for the designer to obtain only 
statistical information concerning thé response of a communication system to the 
Signal portion of the input. 


Saturation Due to Noise 


It is generally quite important that a noise component of the reference 
input not produce excessive saturation of any of the stages of the feedback- 
control system. Calculation of the degree of saturation by statistical tech- 
niques is extremely difficult in the general case because it requires a knowl- 
edge of the amplitude distribution of the variable rather than merely its rms 
value. On the other hand, Newton“ shows that one can get a reasonable engineer- 
ing solution to this problem quite readily if one can assume that the noise input 
has a Gaussian distribution. 


If the noise has a Gaussian distribution of amplitudes and its average value 
is zero, the density of the noise amplitudes is given by the well-known 
expression 


ie 
Wix,) = Eb o7 (2/2) (x, /%) (2) 
271 OF 


where x, represents a given amplitude of noise and 0}, the rms value of the noise. 
The probability P that the magnitude of noise does not exceed a given amplitude 
A is thus equal to the area under the density curve, between the limits 


-A < Xy < - As ie€e, 
+A 
P =a dx W(x) (2) 


-A 


If A represents the saturation level of an element through which the noise must 
pass, then the probability P, that the noise will saturate the element, i.e., 
the probability that the magnitude of noise will exceed A, is obtained by 
subtracting Eq. (2) from unity. 


s ms 
se es has if dx W(x) : (3) 


oA 


To apply Eq. (3) it is convenient to nondimensionalize the noise amplitudes and 
saturation level in terms of the rms value of the noise. Define the nondimen- 
sional noise amplitude as u, and the nondimensional saturation level as @: 


| ST rae a) 
a HOR | (5) 


Then the probability that the noise will saturate an element having a saturation 
level equal to @ times the rms value of the noise is 


| 

P = ai =e 

eels du_ eon (6) 
s SS i n 


A plot of this probability is given in Fig. 1, as a function of the nondimen- 
sional saturation level cd. 


For Fig. 1 to apply to any variable in the system, it is necessary that that 
variable have a Gaussian distribution in amplitudes. Generally, one can assume 
this condition holds if he assumes that the input is Gaussian and that the per= 
centage of time the system is in saturation is quite small because (a) if the 


P, * PROBABILITY OF NOISE SATURATION 


aks SATURATION LEVEL 
RMS NOISE 


Fig. 1 - Plot for determining satura- 
tion due to noise. 


per-unit time of saturation is small, it is logical to assume that the system 
behaves essentially as a linear system and (b) if the input to a linear system 
is Gaussian, the system response must be Gaussian. 


Figure 1 can be employed as follows. Assume that the percentage time of 
saturation is to be limited to some value, say 5 per cent. Figure 1 shows that 
for a probability of saturation of 5 per cent the system must pass linearly noise 
amplitudes up to 1.96 times the rms value of the signal. Thus, to determine the 
noise level that must be passed linearly at any point in the system (in order to 
limit the probability of saturation by the noise, at any point, to the given 
value of 5 per cent), the rms value of the noise at that point should be computed 
on a linear basis and mltiplied by the factor 1.96 to obtain the required linear 


range. 


. , idered in 
On the other hand, the signal as well as the noise mst be cons 
determining the required saturation level of the stages in a system. Adding the 
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noise amplitude which must be passed linearly to a time-varying plot of the 
Signal component at a given stage gives a plot of the signal-plus-noise to be 
passed. In any particular case, engineering judgment can be used to determine 
whether the saturation level of that stage must be greater than the maximum value 
of this signal-plus-noise plot or whether a lower value can be allowed. There is 
a fundamental difference between the effects of saturation due to the noise por- 
tion and saturation due to the signal portion. Since the noise is of high fre- 
quency, it drives the system ‘into saturation only for very short instants of 
time; whereas the signal portion, being of low frequency, could maintain satura- 
tion for such long periods of time that the system may not behave in a quasi- 
linear fashion even though the percentage of time during saturation is small. 


BASIC EQUATIONS 


The analysis to be presented is based upon the convolution integral 
described in detail in Appendix A. This integral gives the output time response 
Xo(t) for a given input x,(t) as 


10.8) 
alto if dt h(t) x(t - 7) (7) 
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where h(@) is the impulse response of the system. The impulse response repre- 
sents the inverse Laplace transform of the system transfer function H(s), defined 
as 


X, (s) | 
H(s) = ts) (8) 


where X,(s) and X;(s) are the transforms of x,(t) and x,(t) respectively. 


There are four statistical functions pertaining to the variables x.(t) and 
X(t), which are as follows: 6;;(T), the autocorrelation function of the input 
x4(t); %9(T), the autocorrelation function of the output x,(*%); @,,(F) and 

( i), the cross-correlation functions of the input and the output. The equa- 
tions for these functions are given below. 


+T 
B5,(7) = sae o- i dt x, (t) x, (t 28) (9) 
— -T 
+T 
Ppcelsr vie ee a fis dt x(t) x(t + 1) (10) 
: -T : : 
G5 c(t) = tard = 3, if dt x, (t) x(t + 1) 3 (11) 
? ray i t 


ue 


647) = scene a i dt x(t) x, (t + 7) (12) 
-T 


Techniques for computing these functions from input and output data are described 
in Appendix A. 


Some important characteristics of the correlation functions can be obtained 
by inspection. The equations for the autocorrelation functions 034) and 09,(T) 
have the same value if t is replaced by -7; so that 


9..(t) = $..(-7) (13) 
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Thus, the autocorrelation functions are even functions of T; the cross-correla~ 
tion functions, on the other hand, are generally not. If in Eq. (11) the varia- 
ble T in 0;,(T).is replaced by -7, the relation becomes equal to Eq. (12) for the 
other cross-correlation function 0,5;(7). Thus, 


Pomme Pt 0) (15) 


If the autocorrelation function 03, (7) in Eq. (9) is calculated for t= 0, it 
becomes 


+T 
oe . ae 2 
PogtOhe pre aot f dt x," (t) (16) 
— 00 -T 


This is by definition the average of x; ? orn other words, the mean-square 
value of the input, which is designated as x;°. Thus, 


0,5 (0) = x52 (17) 
and, likewise, 


¢ el Oy)sce. a & (18) 


Oo 


EXPRESSION RELATING MEAN-SQUARE VALUE OF 
OUTPUT TO STATISTICAL CHARACTERISTICS OF INPUT 


As has been shown, it is desirable to be able to calculate the mean-square 
value of the output response of a system from a statistical condensation of the 
input, rather than from the complete input plot. To find how this calculation 
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may be performed, consider the definition of the uote Corbet value of the output, 
which is 


Sie. , uf 
a eee A ae if dt x(t) x(t) (19) 
T —® 00 ~T 


Substitute into Eq. (19) for each of the variables x,(t) the convolution integral 
of Eq. (7), using separate time variables T7 and 72 6 differentiate between the 
two X,(t) variables. 


+T oo A. 
ne ; s > —_ 
Xoo = shoe oT fe fer, h(T}) x, (t-T, ) f dT, h(To) x4(t ia) (20) 


Reverse the order of integration. 


oS eo +T 
ae = f at n(T,) [ at h(T 5) ss or sh dt x, (t-7, ) a(t) | (21) 
-T 


ie) Oo 


The expression within the brackets may be put in the form 


+T 
Lim oa f dt x, (t ae AS [te ~ T,) a Poe, - to) (22) 
T —®e oo “? 


As shown by Eq. (9), this is the autocorrelation function of the input with 
(T] -T2) representing the time shift; i.e., it is 0;,;(€] - To). Thus, Eq. (21) 
may be expressed as 


00 
ae Se f da, h(t, ) dt, h(t, ) ie AL =T >) (23) 
, 


oY 9 


Equation (23) expresses the mean-square value of the output response of a 
system as a function of the autocorrelation function of the input and the system 
impulse response. Thus, the equation shows that the autocorrelation function of 
a signal contains all the information concerning that signal that is needed for 
determining the mean-square value of the response to that signal of any linear 
system. Although Eq. (23) might appear rather complex to the reader, it can be 
solved readily by transient and transform techniques. One can show from the 
equation that the mean-square value of the output can be found by treating the 
input autocorrelation functions as a transient input to the system. On the other 
hand, it is easier to demonstrate this by directly treating the input autocorre- 
lation function as a transient input and calculating the response of the system 
to that transient input. 
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RELATING CORRELATION FUNCTIONS BY TRANSIENT APPROACH 


The correlation functions are functions of the time shift T but can be 
conveniently considered also as functions of real time t, and the output auto- 
correlation function 0 o9(t) can be calculated by transient techniques from the 


auput autocorrelation function 0;;(t). The mean-square value of the output 
X° can then be found by determining the value of 0,,(t) at zeroT. 


Consider 05;(t,) as a function of the real-time variable 7) and an input to 


a system with the impulse response h(t). The time-response of the system to this 
transient input is given by the convolution integral of Eq. (7) as 


[e,@) 
if dt, h(t,) #,,(7, - T,) (2h) 
(0) 


Using the general form for the input autocorrelation function given in Eq. (9), 
substitute into this integral the expression for $,;,(] - To). 


co 


+ 
iL 
of af R(t) Tim 3 if at x, (t) x, [t (CHE T)] (25) 
Reverse the order of integration. 
+T ae 

1 = 26) 
bin et ue dt x; (t) Sar; h(T>) x, |(t +7) -T2 ( 
T—» co Adak -T = , J 2 1 [ 1 | 


_ The expression within the large brackets is a convolution of the input x3 (t+T,) 
with the system impulse response ho) and, hence, represents the output 
X(t +7), since by Eq. Cie 


00 
x i (t + T,) = i, dt, h(t, ) x, {(t Ge me = | (27) 


re) 


Substituting Eq. (27) into Eq. (26) gives 


or Ah 
; 1 
Lim oT iE dt x, (t) x(t + 14) (28) 
T —e 00 _T 


oar: . lation 
ring this with Eq. (11) shows that it represents the cross-corre 
i eee rN (T)- Thue, if the input autocorrelation function 055 (7) is taken 
asa transient time-input to the system, the system time-response to this input 


is the cross-correlation function 0; 9(7). 
) ; %) about the zero time-shift 
rotating the cross-correlation function 9; ¢ | 
axis "The cross-correlation function ),;(t) is formed. Now determine the system 
pesponse if 05;(¢) is taken as a transient input to the system. The convolution 
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integral of Eq. (7) gives as tne system time-response to the transient input 

O04 (T7)5 
(0, 8) 
f Ce 9 a ue et (29) 
18) 


Substitute into Eq. (29) the expression for 0); (Tj - 73) obtained from the 
general definition of Eq. (12). : 


foe) met 

‘ AL: 
fan. h(t.) ae oT / dt x(t) x, |t + (tT, ~ ) (30) 
ro) ra ~T 


Reverse the order of integration. 


+T 
Lim oad fae x(t) LJ dt, h(t.) x [tt Ait so -3]| (31) 


T—# 00 
-T fe) 


co 


The expression within the large brackets is the convolution of the input 
x4 (t + ©) with the system impulse response h(T3) and hence represents the 
output x(t + 71). Thus, Eq. (31) becomes 


+7 
j 1 
Lim oT i dt x(t) x(t a T,) (32) 
T —® oo _T 


which is the output autocorrelation function O50 (Cq)> given in Eq. (10). Thus, 
if the cross-correlation function 53 (F) is taken as a transient time-input to 
the system, the system-time response to this input is the output autocorrelation 
function )9,(7). 


The complete transient method for calculating 0,,(7) from 0;;(%) is illus- 
trated in Fig. 2. The input autocorrelation function 055) is be in Fig. 
2(a). It is symmetric about zero T, and its value at 7 = 0 is the mean-square 
value of the input x;2. With 03;;(7) as the transient input to the system H(s), 
the output is 0;,(7), which in general is not symmetric about zero 7. 


Rotating the plot of 0;,(7) about the zero 7 axis gives the plot of 0); (T) 
shown in Fig. 2(b). With 0,;(7) as the transient input to the system H(s), the 
output is YoolZ), which is symmetric about the zero t axis. The value of 0,,(t) 
at zero T is x5*, the mean-square value of the output. 


The transient calculations described above for obtaining O50 (T) from 054) 
can be performed readily by the graphical method.+ It is also possible to per- 
form them by an analog computer, but this may be quite difficult because it 
requires that input-table plots of both 944(€) and 09;(T) be prepared. Conse- 
quently, when an analog computer is to be employed, it is often best to modify 
the calculations as described in Appendix C so that only an input-table plot of 
0345(€) need be prepared. | 
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,, (vr) an (z) 


i'm ae 
INPUT OUTPUT 


(b) RESPONSE TO ¢,( 


Fig. 2 - Transient response to correlation 
functions. 


RELATING CORRELATION FUNCTIONS BY TRANSFORM METHOD 


The mean-square value of the output may be calculated by transform tech- 
niques. In this section, the following basic equation is developed for the mean- 
Square value of the output. 


ro) 
a es 
pa OF da 


where H(jw) is the transfer function of the system and ¥;;(@) is the Fourier 
transform of the input autocorrelation function 0..; 4 (7). Thus, in the transform 
procedure for calculating the mean-square value of the output, the transform of 
the input autocorrelation function is calculated and multiplied by the square of 
the magnitude plot of the system transfer function; the mean-square output is 
then proportional to the area under the resultant curve. 


H(ja) | 044.) (33) 


Empirically computed autocorrelation functions are necessarily of finite 
length, and the finite length is an important cause of error in the transform 
method. Techniques for Pee ome ae transform computation in a way to minimize 
this error are given by Ross.- The finite length of the autocorrelation function 
also causes errors in the transient method, but it appears that the errors should 
be less severe when the transient method is used. However, a detailed study of 
the problem is required before this conclusion can be substantiated adequately. 


ae: 


The main advantage of the transform approach over the transient approach is 
that the transform of a correlation function expresses the data in terms of its 
frequency information, which often is more convenient to handle in an engineering 
analysis. 


Since the correlation functions are not zero for negative values of T, it is 
necessary to transform them by the Fourier transform rather than by the Laplace 
transform. The Fourier transform and its inverse are as follows. 


+ 00 
Bue) = J Eifion awe dt (3h) 
~00 
+00 
te - jot 
£(t)-= 27 F(w) e 9 dw (35) 
~00 


The factor (27) is sometimes placed with the F(W) integral rather than with the 
f(t) integral, but this makes no difference in the over-all calculations provided 
the pair of equations is consistent. 


The Fourier transforms of the autocorrelation functions 0;4;(T) and 0, &) 
are called the spectral densities or, sometimes, power spectra and are dest enated 
as 03;(W) and 0,,@). Thus, 


054(©) = J 055(&) e7I** at (36) 
05;() = os fest 5" aca, (37) 


The Fourier transforms of the cross-correlation functions 0, and 0,;(T) in 
like manner are designated at 0; (©) and 4 (w), respectively. 


It is possible to calculate the spectral density by transforming the time 
function directly, rather than by transforming its autocorrelation function. 
Designate the Fourier transform of the time function x(t) measured over the 
time interval -T < t< + T as Xp(~). The spectral density of x(t), designated 
0(), can than be shown to be equal to 


O(o) = Lim Fe Xp(w) K_*(o) (38) 
T+00 
stim be |aq¢e |’ 
2 oF Xp fia (39) 
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Equation (39) is proved in Appendix B. Note that the asterisk in Eqe (38) 
designates the conjugate of Xm 


. Theoretically, the spectral densities could be measured directly from the 
input and output signals by means of a narrow bandwidth, variable-frequency fil- 
ter. If the input signal x; is passed through a narrow bandwidth filter centered 
at the frequency ®,, the average value of the output from the filter, divided by 
the bandwidth of the filter, approaches 0;;@) as the bandwidth approaches zero. 
On the other hand, this method of measurement is often quite difficult to perform 
in practice and is subject to error because of nonideal characteristics of the 
filter. However, the concept is important primarily because of the physical 
significance it gives to the spectral density. 


Since the time response of the system to the transient input 0,,(v) is 
Oio(%), the transforms of these correlation functions must be Soe by 


0,0) = H(5) 055) (40) 


where H(jw) is the system frequency response. Similarly, since the time-response 
of the system to the transient input 0,;(T) is 0,,(¢), then 


0,9) = H(36) 0,3). (1) 


Now, it can be shown readily from the basic Fourier transform relations in 
Eqs. (34) and (35) that 


F(-a) = F*@) (2) 
and if fo(t) is equal to f)(-t), then the transforms are related by 
Fo@) = Fi (-4). (3) 
Therefore, the cross-correlation transforms are related to each other by 
M3) = Oo(-4) = Oy) (lui) 
since 09;(t) is equal to ¢;,(-7). Substitute Eq. (4h) into Eq. (41). 
O59) = H(5@) 0%) (45) 
Substitute for @*;,() the conjugate of Eq. )0. 
Ooo) = H(54) HB (50) O%45@) (46) 
Now, since 0;5(T) is equal to 0;4(-?), then by Eqs. (2) and (13) 
0, ©) = 944(-0) = 5) . | (u7) 
Substitute Eq. (47) into Eq. (h6). 
9) = H(j4) H"(ja) %4@) (48) 
= | 1(30)|? 456) (49) 


f the 
Thus, the output spectral density is equal to the product of the square o 
magnitude of the system transfer function times the input spectral density... 
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There are a mumber of convenient relations between the spectral density and 
the autocorrelation function that can be obtained by inspection from Eqs. (36) 
and (37). Setting w equal to zere in Eq. (36) gives 


034 (0) =f 0550) dt. (50) 


Thus, the de value (zero-frequency value) of the spectral density is equal to the 
area under the autocorrelation function. Setting 7 equal to zero in Eq. (37) | 
gives 


x47 = 045(0) =F (440) a0 (52) 
= f o43(e) af (52) 


where f is the frequency in cycles per second. Thus, the area under the 
spectral-density curve (with f as the frequency variable) is equal to the value 
of the autocorrelation function at zero 7, which is the mean-square value of the 
signal. To calculate the mean-square value of the output, apply Eq. (51) to the 
output. 


Xoe F o9(0) = FE (Og) do . (53) 


Substitute for 0),.(@) the expression im Eq. (8). 


52 = Saf |B) : 0;;() do (33) 
=f H( 527) : 0, ,(£) af (54) 


Equation (33) was the basic equation given above for.calculating the mean=square 
output. 


In computing the transforms, it is much easier to express the equations as 
integrations of real variables, rather than as integrations of complex variables. 
This can be performed quite readily with the autocorrelation functions, since 
they are even functions. The transform integral of Eq. (36) can be expressed as 


3; @) = f © (cos WT = j sin wt) dt (55) 
because 
=-jO . ; 
e = cos @=- j sine. : (56) 
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Since 055 (7) is an even function of 7, the sine ee fi ; 
the equation reduces to a P n of Eq. (55) is zero, and 


oo 


05; ©) =f 05; (7) cos wl dt (S7) 
Eis fae cos WT dt. (58) 


{@) 


This integral is always real, and so the phase of the spectral density mst 


eiweye be zero. Besides 0;;() must be symmetric about the zero-frequency axis. 
US 


0;;@) = O14 (-0) = O5,@) . (59) 


Since the spectral density is an even function of frequency, it can be shown 
readily that the inverse transform integral of Eq. (37) reduces to 


0547) = ro) cos wT dw (60) 
pee 
=F fre cos wl dw. (61) 
ae 
APPENDIX A 


Graphical Interpretation of Convolution and Correlation Integrals 
Convolution Integral 


The convolution integral can be derived by physical reasoning. Consider the 
input x;(t) as representing the sum of a large number of pulses of width At and 
magnitude x4(t). If the time interval At is sufficiently small with respect to 
the response time of the system, each pulse has the effect of an impulse of the 
same area, pees Figure 3(a) shows the response of a system to a unit 
impulse. Figure 3(b) shows an input x,(t) to the system and one of the pulses 
into which x;(t) is considered to be divided. The response to this pulse is 
obtained by multiplying the unit impulse response h(T) by [x;(t) At] to obtain 
the response curve shown in Fig. 3(b). At the time t, the response to this pulse 
has the value / 


[24 (t) At] n(ty - t) (62) 


as shown in the figure. The total response x,(t,) to the input x, (t) at the time 
ty is obtained by adding together the responses te all the pulses into which 
x;(t) is divided 


lt) = Da ad) nce - ; (63) 
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0 Te 


(a) UNIT- IMPULSE RESPONSE 


" [xj At}hcr,-t) 


PRODUCT OF THE 
TWO CURVES 


UNIT IMPULSE RESPONSE TIMES 
VALUE OF IMPULSE X;(t) At 


NET AREA= X,(t,) 


(b) CALCULATION OF RESPONSE TO INPUT 


Fig. 3 - Derivation of convolution Fig. h - Alternate interpolation of convo- 
integral. lution process. 


for all the pulses between t = -coand t = tj. Passing to the limit gives 
al 
x(t;) = i dt x, (t) h(t, - t). | (61) 


The convolution process is illustrated in Fig. h in a slightly different 
manner than in Fig. 3. The time tj at which the output is being calculated can 
be considered to represent the time of observation. The unit impulse response 
h(t) is shown projected from the time of observation t] into the direction of 
past (or negative) time and hence represents h(t, - t) as shown. Multiplying 
the input curve x;(t) point by point by the impulse response h(t j-t) gives the 
product curve x;(t) h(t j-t). By Eq. (64) the response x, at present time t, is 
the integral of this product curve from -- up to the present time t, and hence 
represents the net area under the product curve. 


The construction of Fig. shows why the unit impulse response h(7) is 
termed a "weighting function." The impulse response h(t) mltiplies or 
"weights" the past values of the input. For systems with a finite bandwidth 
the impulse response h(t) goes to zero at infinite 7, and hence the values of 
the input x, (t) far in the past must have little effect upon the present output 
Xo(t,). The more recent values of the input are, in general, weighted much more 
than long-past values. It is obvious also that values of the input which have 
not yet occurred cannot affect the present value of the output. This condition 
is assured in the convolution integral because the impulse response is- zero for 
negative T. 


The variable T described above is equal to 
eo (t) ” t) (65) 
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and therefore can be considered to represent elapsed time, a time variable from 
present time t, projected into the past as shown in Fig. lh}. To express the 
convolution integral in terms of T, note that by Eq. (65) 


dT= —- dt. (66) 
Substituting Eqs. (65) and (66) into Eq. (6h) gives 

(t,-T) = +1 
(4) -5 i (-a7) x, (t,-7) n(7) . (67) 

(t,-T) = - 0c 

Expressing the limits in terms of T gives 
x, (tz) = far h(t) x;(t; -T) . (68) 

fs) 


Equation (68) is the more convenient form of the convolution integral which is 
used in this paper. 


Autocorrelation Function 


The process of computing an autocorrelation function is illustrated in Fig. 
5. Figure 5(a) shows a plot of an input function x.(t) and shows the same plot 
shifted to the right by the time shift Ty to obtain x;(t +%) « These two curves, 
x;(t) and x;(t + 7% ,), are multiplied together point by point to obtain the prod- 
uct curve [x,(t) x;(t + 7)], shown in Fig. 5(b). The average value of this 
product curve represents the autocorrelation function 0;;(% |). Thus, 044(7 1) is 


$.,(7,) = Average {x,(t) x(t + +1) f (69) 
t. 
= Lim hts if x, (t) x,(ttr,) a} (70) 
(t,-t,)—e@ ut | 


In practice, of course, the averaging is always performed over a finite period of 
time (to - t;). However, if this period of time (to - t,) is much greater than 
the maximum tine shift T], the point 0,;(T ,) of the autocorrelation curve may be 
considered to be approximated adequately. 


The expression for 0,;(7,) given in Eq. (70) is equivalent to the more 
common form 


+r 
§,, (7) = Saas o if dt x, (t) x, ( t + 7) (71) 
ha. -T 


The averaging process of the product curve described in Fig. 5(b) gives one 
point @;;(& ) on the convolution curve, as shown sy carry a (os json Oy ae obvious 
that if the interval (to - t,) may be considered infinite, it makes no difference 
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in Fig. 5(a) if the second curve is shifted in the positive or negative direction 
by the amount 7]. In either case there are two identical curves shifted from one 
another by the amount 76 Thus, the autocorrelation function 0 (T) mst have 


the same value at -?U that it has at +7; i.e., 0;4;(%) mast be symmetric about the 
zero T axis. 


If the input seh does not have any periodic components, the variations in 
x;(t) are uncorrelated for infinite T shifts. Consequently, as T approaches 


TIME SHIFT 


TIME ——e— 


(a) INPUT CURVES 


NET AREA (t>-t,) D:(T) 
[Xi (t) Xj (t+7)] E ot) Piilt 


MEAN - SQUARE 
VALUE OF X;(t) 


Fi T) 


AVERAGE VALUE SQUARED 
— 2 
(Xj) 


=f, OT; 
TIME SHIFT T 


(c) AUTOCORRELATION FUNCTION. $(r) 


(X, - Xj)" = MEAN - SQUARED VALUE 
OF X;(t) FOR ZERO AVERAGE VALUE 


te) 
TIME SHIFT T 


(J) AUTOCORRELATION FUNCTION FOR SIGNAL WITH 
ZERO AVERAGE VALUE 


Fig. 5 - Description of autocorrelation procedure. 
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infinity the autocorrelation function 0;;(T) approaches a constant value, and 
that constant value merely represents the square of the Zvepage value of the 
input x1. Symbolically, this final value of 044(t) is (x;)*. 


In practice the final valve of 0;;(7) is almost always of no interest. It 
is therefore removed, either by subtracting from the original input curve x3 (t) 
its average value Xj; or by subtracting from the aytocorrelation function 0; (75) 
the square of the average value of the input (X;) « in either case the auto- 
Dar Raaaioaanae actually used has the form of Fig. 5(d), with a zero final 
value. 


This practice of subtracting out the average value of the input is justi- 
fied as follows. In most cases the average value of the signal considered really 
is zero, and it appears in the computation only because a bias is developed by 
the original measurement equipment or in the device computing the autocorrelation 
function. Besides, if there really is a steady component (average x; ) in the 
input signal, the response of the system to the steady component could be calcu- 
lated best by considering the steady component separately from the rest of the 
input. 


It is apparent that this computation of the autocorrelation function could 
be extremely complicated if performed manually. Consequently, the calculation of 
autocorrelation functions is usually performed by a computer, except for some 
very specialized types of input functions in which analytical calculation is 
possible. Very often only the positive half of the autocorrelation function 
(i.e., for positive T) is given. This is because the negative half is merely its 
mirror image and, hence, is not necessary. Nevertheless, one should always con- 
sider the autocorrelation function as having a symmetric negative portion. In 
this sense, it is quite unlike a system impulse response h(t), which mst be zero 
for negative 7. That is why the Fourier transform is applied to 0;,;), whereas 
the Laplace transform can be used with h(t). 


Cross-Correlation Function 
The cross-correlation function 0;,(¥) is obtained in the same manner as the 


autocorrelation function 0;;(t), except that the input x;(%) is multiplied by the 
shifted output x,(t +7) and the average of the resultant product curve caleu- 


lated. Thus, | 
ss. , +7 
659(%) = ian i UE: dt-x, (t) x,(t_+ 4). (72) 


If the output curve x,(t) is fixed and x(t) shifted forward by T to obtain 
x(t + 7), the other cross-correlation function 05; (%) is obtained, which is 


+T 
Os Cipas cite Ge ff dt xq(t) x(t 4 7). ir C73) 
It is obvious that the same value is obtained by shifting the output curve 


forward in time by Tt, or by shifting the input curve backward in time by 7}. 


Hence, 


2 : : Wo (T) FH Oog(-D. (1) 
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On the other hand, the sense of the relative shift between the input and output 
curves (i.e., which one is shifted forward) does make a difference. Therefore, 
the cross-correlation functions are generally not symmetric about the zero axis; 
ieee, O3,5(t) is generally not equal to $4,(-T). 


In calculating the cross-correlation functions, the average values of both 
the input x;(t) and output x,(t) curves are generally subtracted out, so that the 
final values of the cross-correlation functions approach zero at infinite T (pro- 
vided there is no common periodicity in. the two signals). If these average 
values were not subtracted out, the cross-correlation functions would approach at 
infinite T the value (x;)(x,), which is the product of the average values of the 
two signals x: (t) and x(t). If the cross-correlation function approaches a 
periodic waveform at large UT, this indicates a common periodicity in the two 
Signals. 


APPENDIX B 
Relation Between Spectral Density and Transform of Time Function 
The following is a demonstration of the relation given in Eq. (39), which 
relates the spectral density to the transform of the actual time-varying signal 
rather than to the autocorrelation of that time signal. 
Define the function x(t) as equal to the time signal x(t) over the time 


region -T < t< +T and equal to zero outside that region. The transform of xn (t) 


oO 
Sn radi elie he AR) pe 
— 00 
A correlation function On (7) may be defined in terms of Xm (+) as follows. 
[6 @) 
atk ot Sits [es ory (eh aes ter a at (76) 
-0O 


~The autocorrelation function 0(T) of the total signal x(t) is the limit of Op (T) 
as the time region becomes infinite. 


O(+) = Lim ¢.(7) 


Ai (77) 
The transform of p(T) is designated as 0p(w) and is equal to 
Ono) = f aq Oy(ey) owFA (78) 


Therefore, the spectral density of the complete input x(t), which is designated 
b(), is the limit of 0 (©) as the time region becomes infinite. 


00) = Lim (0) (79) 
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Substitute into Eq. (78) the expression in Eq. (76) for Op (77). 


od 


-3 1 
On) =f e Joe oF fo Xm (t) Xp (t + Tt) (80) 
Now, ea can be expressed as 
Suite eroeee tg) cet (81) 
so that Eq. (80) is equal to 
sat Z's r 
0,() = ~: (Seif er (toe Jo(t + Leber (82) 
Reverse the order of integration. 
j ~jo(t 
0 (©) = = fa EG Xp (t) fam e Jo( DV xn(t4%) (83) 


In terms of the second integral the variable t is a constant; therefore, for that 
integral dt) is equal to 


aT, = a(t +7) . (8h) 


Substitute for (t +71) in the second integral of Eq. (83) the variable 7. 
Equation (83) becomes 


Op(o) = Fe if at 0% xa(t) fo aty oo? x(t) . (85) 


Comparing the two integrals with the expression in Eq. (75) for Xn) shows that 
the second integral is the transform of x(t) which is Xp@), and the first is 
the conjugate of this transform which is Xin() « Thus , 


| Eales 
bp) = L x(0) Xq(0) . (86) 
Since the limit of 0p(@) is the spectral density %(), then O(w) is 


Oo) = Lim 22 Xp) Xp) (87) 


T— 


1 2 
= Lin an) | (88) 


which was to be shown. 
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APPENDIX C 
Method for Computing Output Autocorrelation Functien on an Analog Computer 


In calculating the mean-square value of the output signal for a given input 
autocorrelation function, analog computers can be very useful, especially if a 
great many settings of the system parameters are to be tried. This calculation 
could be performed by an analog computer with a two-step precedure: first, 
determine the response 049(T) of a simulated system to the input 0;;(T); and 
second, reverse the plot of 0;,(€) in time to get 003 (©) and find its response, 
which is 0,,(). However, this procedure may be quite inconvenient because it 
requires that an input-table plot of 053 (8) be prepared for each parameter 
setting of the system. On the other hand, the procedure can be modified as 
described below, so that the output autocorrelation function can be calculated 
with only a plot of 0;;(t) required for the input table. 


Since the. output autocorrelation function 056 (7) is the time response of 


the system with the cross-correlation function 09;(T) as a transient input, then 
0,.(t) can be expressed in terms of 0,;(€) by the convolution integral. 


oD 


oo(ty) = f at nlp) Yoga - Te) (89) 
Since 09;(T) is equal to pecuiciee 
041 - Fo) = bio(Co - %)- | (90) 
Substituting Eqe (90) into Eq. (89) gives 


V0 (71) = [a h(o) 0, (2 - Tt). (91) 
For 71 equal to zero the output aagocorrelati on function 0,,(7;) represents the 
mean-square value of the output x,¢, and Eq. (91) becomes 
ez | 
2 
X= (0) 2 far h(T>) 0455). (92) 
0 


Figure 6 illustrates graphically how Eq. (92) may be computed. The response 
of the system to the transient input 0;;(7) is 0;,(%), shown in Fig. 6(b). The 
response of the system to a unit impulse at T = é is the unit impulse response 
h(@), shown in Fig. 6(c). Multiplying point by point the plots of 05,(T) and 
h(t) gives the curve in Fig. 6(d). The area under this curve is the mean-square 
output xe ° 


Figure 7 shows the computer set-up for computing the mean-square output 
value. On the input table there is a plot of the input autocorrelation function 
054 (T) and a trigger mark placed at the point of zero TU. The input autocorrela- 
tion function 0;;(T) is fed to an analog of the system being studied, and the 
output is bi o(t)- The trigger actuates a pulse generator feeding a sharp pulse 
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Fig. 6 (left) - Steps in 
computing mean-square 
output by analog com- 
puter. 


Pig(T2) 


Fig. 7 (below) .- Analog 
computer set-up for 
computing mean-square 
value. 
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TIME T—~ 
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SYSTEM 
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which is effectively a unit impulse into a second analog of the system, and the 
output is the impulse response h(T). The outputs from the two analog systems 
are fed to a multiplier unit to form the product 6;,(T) h(t) and this is fed 
into an integrator, which delivers in the steady state the mean-square output 
Xy e 


One practical difficulty that the technique has is that it is difficult to 
use impulses in an analog computer without producing computer saturation. A 
simple way to alleviate the problem is to simulate in the impulse-response path 
the transfer function of the system multiplied by the factor s and feed into 
this, in place of the impulse, a unit step actuated by the trigger. 


The mean-square value of the output, calculated from the computer set-up, 
represents the value of the output autocorrelation function for zero Pee 
obtain other points on the output autocorrelation function plot, shift the 
trigger with respect to the input autocorrelation function by the desired time 
shift 7. That this is so can be seen by expressing Eq. (91) in the form 


[—) 


0 Cz) = { at, h(T, +1) 03,2) « (93) 


fe) 


: wane 5 


Thus, a shift of the trigger in the positive direction with respect to 0;;(T) by 
the amount 7; produces the point on the output autocorrelation function at the 
positive time shift Tj. 
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A SURVEY OF TECHNIQUES 
FOR THE ANALYSIS OF SAMPLED-DATA CONTROL SYSTEMS 


Gordon J. Murphy | 
University of Minnesota 


% Ralph D. Ormsby 
Minneapolis -Honeywell Regulator Company 


Summary 


The present use in control systems of pulsed-data links, track-while- 
scan radar, digital computers, and many other intermittently operative devices 
has stimulated interest in the analysis of sampled-data control systems. The 
tremendous effort now being expended to develop the techniques for analyzing 
such systems has resulted in the proposal of several methods of analysis. This 
paper presents an up-to-date discussion of the state of the art and lists pro- 
posed techniques of analysis with an explanation of the applications and limi- 
tations of each. 


Introduction 


The discussion which follows is concerned chiefly with the three princi- 
pal methods of analysis currently exploited in the literature: impulse-response 
analysis, frequency-response analysis, and transfer-function analysis with the 
aid of the z-transformation. 


In the very simple sampled-data feedback control system in Figure 1, the 
error E between desired output R and actual output C is sampled at dis- 
crete moments of time. The sampled error signal is represented by E*. The 
sampling may be either periodic or aperiodic; in the latter case, it may be de- 


Fig. 1 = A simple sampled-data system. 


pendent upon some characteristic of a signal, usually the system error. The 
impulse-response analysis is the only one of the three principal methods which 
is directly applicable to the aperiodically sampling system. 


Impulse Response 


The impulse-response analysis is the most informative but usually the 
most laborious to apply of the three principal methods. !7 Between sampling 
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instants the system in Figure 1 operates as an open-loop control system with 
pulse inputs uniformly spaced intime. The instantaneous output can be approxi- 
mated by a superposition of appropriately weighted impulse responses of the 
open-loop system, each impulse response commencing at the sampling instant 
for which its weighting factor was determined by the error sample. 


The application of this procedure to a simple system with a unit-step 
input is illustrated in Figure 2. The sampling in this system is periodic with 
period T. The impulse at t = 0 is of unit amplitude, the magnitude of the im- 
pulse at the second sampling instant is equal to unity minus the magnitude of 
the response of the open-loop system to a unit impulse T seconds after the 
application of this impulse, etc. 


Frequency Response 


Impulse-response analysis becomes very laborious when the open- 
loop transfer function is of high order. For such systems, if the sampling 
is periodic, an extension of the conventional frequency-response-analysis 


Fig. 2 - Response of a particular sampled-data system to a unit 
step input. 


IMPULSE 
TRAIN 


IMPULSE 


MODULATOR 


Fig. 3 = Representation of the sampling device as an 
inpulse modulator. 


techniques may yield.the desired information with less effort than the impulse- 
response analysis. The analysis in terms of frequency response, based upon 
the concept that the sampling device is an impulse modulator, is illustrated 

in Figure 3, which is a redrawing of Figure 1. 
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When considered in this light, the sampling device is seen to generate 
signals complementary to the complex-frequency components of its input signal. 
For an input at a single frequency, the output of the sampling device is found 
to possess a frequency spectrum containing the input frequency and an infinite 
number of so-called complementary frequencies separated from one another and 
from the input frequency by integral multiples of the sampling frequency. In 
general, a given component of the modulator input and all modulator output sig- 
nals complementary to this component have the same amplitude at the modula- 
tor output. The expression which relates the sampled error to the continuous 
error is . 


+ 00 
x (s+jnn) where = on is the sampling frequency. 


oS 
E T n=-oa@ 


In a closed-loop system the complementary signals in the modulator out- 
put signal are transmitted to the input to the modulator. Because of the filter- 
ing action of the system, the components of the frequency spectrum of the signal 
at the input to the sampling device are not, in general, all of the same amplitude. 
Thus, it is possible to sustain a system output of the desired nature (i.e., one 
with a nonperiodic Fourier transform) even though the Fourier transform of 
the output of the sampling device must be periodic. 


A complete frequency-response analysis of such a system must, of course, 
take into account the infinite number of complementary signals; as a result, an 
infinite series is involved. For practical purposes, however, use of a finite 
number of terms to approximate the infinite series yields satisfactory results. 
For most frequency-response analyses, in which extreme accuracy is not re- 
quired, all that need be considered are the system input signal and the comple- 
mentary signal at a frequency equal to the frequency of the input signal minus 
the sampling frequency, since the magnitude of the transfer function at the other 
complementary frequencies is relatively small. If the sampling frequency is 
high and the filtering is fairly effective, satisfactory results can be obtained by 
ignoring all the complementary signals. 


The system used to illustrate the technique of impulse-response analysis 
is used again to illustrate the technique of frequency-response analysis. Curvel 
of Figure 4 is the Nyquist diagram for G, with a change in scale factor to com- 
pensate for the action of the sampling device. By a simple process of vector 
addition, this curve is modified graphically to account for a finite number of | 
complementary signals. Because of the periodicity and the symmetry of the 
function being approximated, the frequency-response curves need be plotted 
only over a range of frequencies equal to half of the sampling frequency. Only 
the first two terms of the infinite series are plotted in curve 2 of Figure 4. 
Curve 3, based upon the first two terms of the series, represents the same 
system with a lower sampling rate. | 


The Bode Diagrams corresponding to curves 2 and 3 in Figure 4 are 
curves 1 and 2, respectively, in Figure 5; and the corresponding Nichols 
curves are curves 1 and 2, respectively, in Figure 6. It is worth noting that 
the question of absolute stability can be answered by application of the con- 
ventional rules to the curves in Figures 4, 5, and6. It should also be noted 


that constant-M contours lose their significance when sampling is involved. 


: 
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Fig.e h = Nyquist diagrams for the sys- 
tem of figure 2. 
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Fig. 5 = Bode diagrams for the system of 
figure as 
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Fig. 6 - Nichols charts for the system of 
figure 2. 


However, a limited amount of information concerning relative stability is con- 
tained in the nearness of the frequency-response curve to the critical point, 
i.e., (-1, 0) in the KG-plane. 


The Z Transformation 


Frequency-response techniques lead to a rapid approximate picture of 
the behavior of the sampled-data system. Whenmore accurate information 
is required and the system is too complex for convenient use of the impulse- 
response analysis, it is advisable to write the exact transfer function for the 
sampled-data system in closed form and apply the techniques of transfer - 
function analysis. Because of the sampling, the transfer function is not a 
rational function of the Laplace transform variable s; asaresult, the use of 
conventional techniques in the s-plane is ruled out. It is possible, however, 
to introduce a transformation, 10 replacing esT with z, where T is the period 
of the sampling device. This transformation results in a transfer function 
which is a rational function of z. Consequently, analysis on the basis of pole- 
zero configurations in the z-plane for sampled-data systems is similar to 
analysis on the basis of pole-zero configurations in the s-plane for systems 
which do not contain sampling devices. 


The z-transform analysis procedure suffers from two limitations: (1)a 
new concept must be mastered by the engineer before proceeding with the 
system analysis: (2) without modification, the z-transform approach provides 
information about the signals at the sampling instants only.as does the frequency 
response method. It is possible to modify the basic method to obtain the value 
of a given signal at other than the sampling moments, but this modification 
requires considerable additional work. 


The weight placed upon these limitations depends primarily upon the 
application. The engineer who spends much time working with sampled-data 
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systems will find it greatly to his advantage to master the convenient and cir- 
cumspective z-transform techniques. The lack of information between sampling 
instants is not normally a serious handicap; the signal is usually sufficiently 
smooth to insure satisfactory behavior between sampling instants when the re- 
sponse at the sampling instants is satisfactory. The following example is 
included to illustrate the method of obtaining the z-transform of a known function 
of time. 


For the function 
f(t) = e 


with sampling assumed to begin at t = 0, the magnitude of the first sample 

is unity; that of the second sample is e~@T; that of the third sample is e-2aT; 
etc. By associating the magnitude of the sample at a given instant with the 
magnitude of an impulse assumed to occur at that instant, it is possible to write: 


a 2aT 


ae -aT ; mae - 
f (t) = Uy (t) + e u(t T)+e uy, {t 25) ; 


The Laplace transform of this signal is 
* 2% - = - - 
Lf (t)=F (ite pies eaaechien so tel 4a aceene 


1 
1 - age oars 


—_ 
= 


Replacing eels with z then gives a function which is said to be the z-trans- 
form of f(t) and which for simplicity of notation is denoted with F(z). Thus, 


F(z) = 
RSS aT 


Note that the z-transform is the same for all functions of time which are equal 
to one another at each sampling instant, even though these functions may differ 
greatly from one another between sampling moments. 


The transfer function of a sampled-data device is defined to be the ratio 
of the transform of the sampled output to the transform of the corresponding 
sampled input. For the system in Figure 7, if the input r isa pulse at 
t = 0 of unit height, with duration less than T, the outputc is e-at. The | 
z-transform of e~@t has been shown above to be z/(z - e-aT), and the z-transform 


of the unit-pulse input is obviously unity. Hence, the transfer function of the: 
system in Figure 7 is 


Fig. 7 - Illustration of the significance of 
the Z-transformed system function. 
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The stability analysis of a closed-loop sampled-data system such as that in 
Figure 1 may be carried out in terms of the Nyquist diagram for C(z)/E(z) 

by plotting this complex quantity for values of z = eIWIT, (Such a curve is ac- 
tually the exact curve which was approximated with a finite number of terms in 
the so-called frequency-response method discussed previously.) However, as 
in the study of continuous-data systems,much more information about the per- 
formance of the closed-loop sampled-data system can be obtained by plotting 
root loci in the z-plane and correlating the pole-zero constellation of. the closed- 
loop system with its transient response. 


The choice of compensation for sampled-data control systems can be 
made on the basis of the reshaping of the root locus in the z-plane. The use of 
a digital computer as a compensating device results in more flexibility in com- 
pensation than is possible with passive networks and permits the design pro- 
cedure to be carried out completely in the z-plane. 


The block diagram of a simple system which includes a digital com- 
puter is shown in Figure 8. The closed-loop transfer function is given by 
the expression 


c S D(z) G(z) 
Re Selene 


DIGITAL COMPUTER 


Fig. 8 = A simple feedback control system which incorporates a 
digital computer. 


The fact that the digital program and the continuous part of the system 
are isolated by samplers permits the multiplication of the individual trans - 
fer functions in the z-domain. The ease with which compensation of con- 
tinuous systems may be designed in the s-plane stems from the ability to 
multiply the transfer functions (the poles and zeros of the combined transfer 
function comprise all of the poles and zeros, respectively, of the individual 
transfer functions). 


The design of digital programs for the shaping of the root locus in the 
z-plane is actually simpler than the design of compensation networks for ane 
tinuous systems. Compensation for continuous systems is commonly obtaine 
from networks composed of resistors and capacitors; the realization of such 
networks is restricted by the requirements that poles be simple and lie on 
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the negative real axis, and that the gain be limited by the response charac~- 
teristics desired. The only requirement for realizability of a digital program 
is that it be expressible as a rational algebraic function of the delay operator z, 
and that b, of the expression below be other than zero. The transfer function 
is of the form 


a +a,(1/z)+a (1/27) +a (1/2) +. sta (1/2) 


D(z) = Rone aC RETIN OME a cee 
bo +b, (i/z) +b, (1/24) + by (1/B) +. . + by (1/zA) 
where the a's and b's are real coefficients. 
The transfer function D(z) may be expressed as 0(z)/I(z) where 0(z) 


is the computer output and I(z) the computer input. The equation above may 
be written as 


LL 2 Oot ak ales ehe) pit 
a) eae ees 2.05 Pg say - +b, ~ u 1 faye Z 
Z Z Z 
1 1 1 
TB ot 5, ag 3 + + an “i 
Z Z Z 


Choosing the a's and b's so that by becomes unity and rearranging give 


s ioe Cres oi el bee ale 
O(z) = I(z) [es 4 Shs + ay 5 + , aes m= 0 (z) |b} = 


a Zz 
1 1 
+bi ——+ ..+b! 
2 2 a 


The power of z distinguishes the sample, and the a's and b's represent 
the weighting of the sample. For example, I(z) a, (1/z2) specifies that the 
input two sampling periods earlier should be multiplied by a,. If bo had 
been zero, the program would have called for a future value ‘of input and 
thus become unrealizable. 


Discussion 


Of the three most prominent methods of analyzing sampled-data control 
systems, the impulse-response analysis is the most straightforward. More- 
over, it provides the maximum information concerning the system response to 
an arbitrary input. Because the application of this method is laborious, it is 
not practical for analysis of the more complex systems. 


A method of analysis which is much simpler to apply but which yields 
less information is the so-called frequency-response method, which consists © 
of a graphical approximation of the Nyquist diagram for the sampled-data 
system. When the only information required pertains to absolute stability or 
rough measurement of relative stability, the frequency-responee analysis is 
rapid and adequate. It is also useful in design work by giving an indication of 
the type of compensation required. For a detailed analysis uf the more com- 
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plex systems, however, transfer-function analysis in terms of the z-transform 
operator is the best of the three methods; it provides more information than the 
frequency-response analysis and is less laborious to apply than the impulse- 
response analysis. 


The three methods discussed are not the only methods which have been pro- 
posed, but they appear to be the most widely accepted. A familiarity with the 
details of these three techniques should enable one to read, understand, and 


apply the results of forthcoming papers on the subject of sampled-data feedback- 
control systems. 
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IRE Standards on Graphical and Letter Symbols 
for Feedback Control Systems, 1955* 


I. INTRODUCTION 


1.1 Standardization of symbols is considered important 
for the exposition of feedback control concepts. The 
purpose is to establish forms for representing letter sym- 
bols and graphical symbols used in block diagrams. 

In the preparation of this standard, it was found that 
there was no complete and self-consistent set of symbols 
in use that appeared to fulfill the requirements of the 
IRE. It was also found that there is a wide variation 
in the symbols used by different industries and pro- 
fessional societies so that it is difficult to choose symbols 
with a universal acceptance in all fields. Therefore, a 
compromise has been made among existing symbols 
used in the electrical field. 


II. GRAPHICAL SYMBOLS FOR BLOCK DIAGRAMS 
2.1 Transfer Element 


2.1.1 A transfer element represents the functional 
relationships (gi2) between a single input signal 
(x1) and a single output signal (x2), in which the 
input signal, indicated by the arrow, is the inde- 
pendent variable. 


2.2 Mixing Point 


2.2.1 The indicated relationship is x3=f(x1, x2). 


2.3 Summing Point 


2.3.1 The indicated relationships are x;= x, —x2 and 


Xg=X1—Xe-+x3—x4. A summing point is a special 


case of the mixing point and indicates the algebraic 
addition of two or more signals to produce one out- 
put signal. An algebraic sign should be indicated at 
the arrowhead for each signal to be added. If the 
number of input signals to be added is large the 
rectangular symbol should be used. 


*Reprinted from the Proceedings of the IRE, 


November, 1955. 
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2.4 Multiplication Point 


2.4.1 The indicated relationship is x;=xyx2. A 
multiplication point is a special case of the mixing 
point. 


2.5 Branch Point 


ia Lata 


2.5.1 A branch point, which indicates that a signal 
is distributed to two or more points in a block dia- 
gram is represented by a heavy dot. Example: 


2.6 Graphical symbols added to a block diagram for 
mathematical purposes shall be shown dotted to 
indicate that they do not represent components of 


the physical system. 


III. STANDARD FOR LETTER SYMBOLS 


3.1 Essential Features of the System of Symbols 


3.1.1 Signals 


Signals are represented by a single letter symbol 
with a single subscript denoting its physical or 
mathematical meaning. The letter x has been 
chosen as the preferred symbol for generalized 
signals. Lower case represents the time domain. 
Upper case represents the complex frequency do- 
main. 


3.2 Transfer Functions 


Transfer functions are represented by a single 
letter symbol with a double subscript, the first 
letter or number of which is the subscript of the 
symbol for the input signal and the second of which 
is the subscript of the symbol for the output signal. 
The symbol g has been chosen. Lower case repre- 
sents the time domain. Upper case represents the 
complex frequency domain. 


IV. EXAMPLES 7 
4.1 Application of the standard graphical symbols and 


the standard form for letter symbols are illustrated in 
the typical block diagrams below: 


IRE Standards on Terminology for 
Feedback Control Systems, 1955* 


1. INTRODUCTION 


1.1. Prior to the preparation of this standard, no com- 
plete and self-consistent set of terminology appeared to 
fulfill the requirements of the IRE. Because of a wide 
variation in the terminology used by different industries 
and professional societies it was difficult to choose a ter- 
minology acceptable to all fields. Therefore, a compro- 
mise has been made among existing terminologies used 
in the electrical field. 


2.. DEFINITIONS 
2.1 Signals 


2.1.1. Loop Input Signal. An external signal applied to 
a feedback control loop. 


2.1.2. Loop Output Signal. The controlled signal ex- 
tracted from a feedback control loop. 


2.1.3. Loop Feedback Signal. The signal derived as a 
function of the loop output signal and fed back to the 
mixing point for control purposes. 


2.1.4. Loop Actuating Signal. The signal derived by 
mixing the loop input signal and the loop feedback 
signal. 


2.1.5. Loop Error Signal. The loop actuating signal in 
those cases in which it is the loop error (See 2.4.10). 


2.1.6. Loop Return Signal. The signal returned via a 
feedback control loop to a summing point, in response 
to a loop input signal applied to that summing point, 
and subtracting from the loop input signal. 

Note: The loop return signal is a specific type of 
loop input signal. 


2.1.7. Loop Difference Signal. The output signal from 
a summing point of a feedback control loop produced 
by a particular loop input signal applied to.that sum- 
ming point. 


*Reprinted from the Proceedings of the IRE, 
January, 1956. ; 


Note: The loop difference signal is a specific type of 
loop actuating signal. 


2.2. Points and Paths 


2.2.1. Mixing Point. In a block diagram of a feedback 
control loop, a symbol indicating the relationship of 
one output to two or more inputs, such that the value 
of the output at any instant is a function of the values. 
of the inputs at that instant. 

Note: If a mixing device in practice contains dynamic 
elements, these shall be considered as transfer elements 
in one or more of the signal paths entering or leaving 
the mixing point. 


2.2.2. Summing Point. A mixing point whose output 
is obtained by addition, with prescribed signs, of its 
inputs. 


2.2.3. Multiplication Point. A mixing point whose out- 
put is obtained by multiplication of its inputs. 


2.2.4. Forward Path. In a feedback control loop, the 
transmission path from the loop actuating signal to the 
loop output signal. 


2.2.5. Feedback Path. In a feedback control loop, the 
transmission path from the loop output signal to the 
loop feedback signal. 


2.2.6. Through Path. In a feedback control loop, 
transmission path from the loop input signal to 
loop output signal. 


the 
the 


2.3. Transfer Functions 


2.3.1. Transfer Function. A relationship between one 
system variable and another that enables the second 
variable to be determined from the first. 


2.3.2. Transfer Ratio. The transfer function from one 
system variable to another in a linear system, expressed 
as the ratio of the Laplace transform of the second 
variable to the Laplace transform of the first variable, 
assuming zero initial conditions. 


( 
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2.3.3. Loop Transfer Function. The transfer function 
of the transmission path formed by opening and proper- 
ly terminating a feedback loop. 

Note: One example of proper termination is a zero 
impedance generator driving the opened loop, and an 
output termination for the opened loop equal to the 
impedance facing the generator. 


2.3.4. Loop Transfer Ratio. The transfer ratio of a loop 
return signal to the corresponding loop difference signal. 


2.3.5. Forward Transfer Function. In a feedback con- 
trol loop, the transfer function of the forward path. 


2.3.6. Feedback Transfer Function. In a feedback 
control loop, the transfer function of the feedback path. 


2.3.7. Return Transfer Function. In a feedback control 
loop, the transfer function which relates a loop return 
signal to the corresponding loop input signal. 


2.3.8. Difference Transfer Function. In a feedback 
control loop, the transfer function which relates a loop 
difference signal to the corresponding loop input signal. 


2.3.9. Through Transfer Function. In a feedback con- 
‘trol loop, the transfer function of the through path. 


2.3.10. Actuating Transfer Function. In a feedback 
control loop, the transfer function which relates a loop 
actuating signal to the corresponding loop input signal. 


2.4. General 


2.4.1. Feedback Control Loop. A closed transmission 
path, which includes an active transducer and which 
consists of a forward path, a feedback path, and one 
or more mixing points arranged to maintain a prescribed 
relationship between the loop input signal and the loop 
output signal. 


2.4.2. Feedback Control System. A control system, 
comprising one or more feedback control loops, which 
combines functions of the controlled signals with func- 
tions of the commands to tend to maintain prescribed 
relationships between the commands and the controlled 
signals. 


2.4.3. Feedback Control System, Linear. A feedback 
control system in which the relationships between the 
pertinent measures of the system signals are linear. 


2.4.4. Feedback Control System, Nonlinear. A feedback 
control system in which the relationships between the 
pertinent measures of the system input and output sig- 
nals cannot be adequately described by linear means. 

Note: A system can be either quasi-linear or nonlinear, 
depending upon operating conditions and performance 
_ requirements. 


2.4.5. Feedback Control System, Quasi-Linear. A feed- 
back control system in which the relationships between 
the pertinent measures of the system input and output 
signals are substantially linear despite the existence of 
nonlinear elements. 

Note: A system can be either quasi-linear or non- 
linear, depending upon operating conditions and per- 
formance requirements. 


2.4.6. Feedback Regulator. A feedback control system 
which tends to maintain a prescribed relationship be- 
tween certain system signals and other pre-determined 
quantities. 

Note 1: This definition is intended to point out the 
fact that some of the system signals in a regulator are 
adjustable reference signals. 

Note 2: It should be noted that servomechanism and 
regulator are not mutually exclusive terms; their applica- 
tion to a particular system will depend on the method of 
operation of that system. 


2.4.7. Servomechanism. A feedback control system in 
which one or more of the system signals represent me- 
chanical motion. 

Note: It should be noted that servomechanism and 
regulator are not mutually exclusive terms; their applica- 
tion to a particular system will depend on the method of 
operation of that system. 


2.4.8. Command. An independent signal from which the 
dependent signals are controlled according to the pre- 
scribed system relationships. 


2.4.9. Disturbance. An undesired command. 


2.4.10. Loop Error. The desired value minus the actual 
value of the loop output signal. 


LOOP 
OUTPUT 
SIGNAL 


LOOP 
INPUT 
SIGNAL 


LOOP 
ACTUATING 
SIGNAL 


FORWARD 
TRANSFER 
FUNCTION 


FORWARD PATH 


LOOP 
FEEDBACK 
SIGNAL 


FEEOBACK 
TRANSFER 
FUNCTION 


FEEDBACK PATH 
Fig. 1 


3. EXAMPLE 


3.1. Application of the standard terminology to a typi- 
cal block diagram is shown in Fig. 1. 
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LETTER TO THE EDITOR 


On Terminology for Feedback Control Systems 


What's all this fuss about terminology for feedback control systems? It is 
well known that when we attempt to design an amplifier an oscillator results. 
Conversely, when an oscillator is desired a highly stable amplifier is realized. 
Why not put this principle to work? A design method for servos is described 
below which makes maximum utility of the above lemma. The approach is dependent 
on the development of a set of inconsistent, contradictory, inherently unstable 
terms to make the servo system believe it is intended for use as an oscillator. 
Such a set of terms is appended to this note. The design procedure then becomes 
a three-step process. 


1. Show the following daffinitions to the systems under design. 


2. Invoke the "Law of the Perversity of Inanimate Objects." We have found 
that this procedure invariably results in a "sailable" product which far 
exceeds our performance spec's. For those people who have little faith in 
this approach we recommend the third step. 


3. Test the system on a psychiatrist's couch (with or without the doubting 
engineer). This simple test can be accomplished without a flight table 
and uses the bare minimum of test equipment. 


In passing, it should be noted that our approach has been extremely success-~ 
ful in designing feedback control loops containing a human operator. The exact 
origin of the attached definitions is unknown. However, when my sister and I are 
released from our present sanctuary we intend to offer a well earned reward to 
the originator of these terms. 


Yours truly, 
Mr. and Miss Informed 
Flying Saucerville 
In Space 
Daffinitions 
Servomechanism: A system, usually of great complexity, in which the output of a 
device is sampled and fed back into the input in order to produce uncontrol- 


lable oscillations. 


Feedback loop: That portion cf a servomechanism circuit which makes possible its | 
instability. 


Follow-up: A device in a servomechanism which is used to follow, as faithfully 
as possible, its oscillations. 


kate-generator: A small but costly electrical device affixed to the output end 
of a servo system for ornamental purposes. 
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Stability: The desired optimum of servo performance, as manifested by violent 
thrashing of an output member. 


Amplifier: An electronic device for superimposing electrical instabilities on 
mechanical. 


Zero back-lash: A whimsical concept much professed by authors of specifications. 


Specification: A form of heroic or epic poetry favored by engineers; a collec-~ 
tion of impossible conditions; a collection of loop-holes loosely held 
together by wistful verbiage. 


Wiring-diagram: A form of the graphic arts which has the characteristic of 
changing form constantly in the manner of cloud formations. 


Nyquist Diagram: An ornamental figure derived from mathematics which demon- 
strates the stability of an obviously unstable system. 


Cathode-ray oscilloscope: A television set devoid of advertizing or wrestling 
matches. 


Production (with respect to servo electronics): The process of painfully assem- 
bling, dismantling, and reassembling electrical components in accordance with 
wiring-diagrams (q.v.)3 a vicious circle. 


Servomechanisms engineer: A soothsayer of mystical powers skilled in mathemati- 
cal rubrics and incantations; a prophet without honor in California or any- 
where else; a dreamer or woolgatherer. 


Performance evaluation: A grateful appreciation of whatever results from a servo 
system. 
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ABSTRACTS OF MINUTES OF THE MEETING OF THE ADMINISTRATIVE COMMITTEE 


OF THE PROFESSIONAL GROUP ON AUTOMATIC CONTROL, INSTITUTE OF RADIO pe one 


HELD AT IRE HEADQUARTERS, SEPTEMBER 17, 1956 


The meeting was called to order at 10:20 a.m. 


PGAG Technical Sessions at IRE 1957 National Convention: Meetings Committee 
Chairman Grabbe has been appointed to represent the Group on the 1957 IRE 
National Convention Technical Program Committee. Mr. David Lindorff will 
serve as his alternate. 


In addition to the papers received at Headquarters in answer to the 
"call for papers" for the Convention, Messrs. Grabbe and Lindorff will 
solicit some invited papers for the PGAC Technical Sessions. It was sug- 
gested that the residue of papers submitted be considered for future use at 
other meetings and that they also be referred to Editor Axelby for considera- 
tion for the Transactions. 


Mr. Grabbe will report to the next meeting of the Committee on the nun- 
ber of technical sessions which will be allotted to the Group by the 
Convention Technical Program Committee and the number and quality of the 
papers received. 


Committee Organization: The chairman reported on his proposed plans for the 
coming year and made the following Committee appointments: 


3a. Awards Committee: Mr. George Biernson was reappointed Chairman of 


this CommLcvee. 


3b. Meetings Committee: Mr. Eugene Grabbe will serve as Chairman of this 
Committee. Mr. David Lindorff will assist him. Mr. Grabbe was 
authorized to appoint other members to the Committee as required to 
procure papers for all meetings. 


3c. Membership Committee: Mr. Azgapetian will continue as Temporary 
Chairman of the Membership Committee until his successor is appointed. 


This Committee will handle all of the following activities: 
membership promotion, chapter promotion, publicity, and student acti- 
vities. Chairman Lozier will appoint four members to handle each of 
the above-named activities, and a new Membership Committee Chairman to 
coordinate these activities. 


3cl. It was moved that the sum of $250.00 be allotted from the 
Group Treasury for a supply of PGAC Membership Brochures. 
Mr. Azgapetian will prepare the text for this brochure and 
submit it to the Chairman and Secretary for approval by Octo- 
ber 17th (unanimous). These brochures will be used for 
membership promotion at the IRE National Convention, WESCON, 
NEC and the Annual Meeting of the Group; for distribution to 
Chapters for their membership activities; and for membership 
ie leh mailings initiated by the Membership Committee.' 
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3d. 


3e. 


3cé. Chapters: Mr. Wilcox agreed to accept the Chairmanship of the 
Wenbersiip Subcommittee on Chapters. He also accepted the 
responsibility for preparing copy for a letter that can be 
used by a Chapter organizer to assist in the formation of a 
Chapter r) 


3c3- Publicity: The primary function of this Subcommittee would be 
to anticipate the necessary publicity for "calls for papers," 
to publicize proposed Chapter activities and interesting fea- 
tures of National and Chapter meetings. This would include 
dissemination of information to the Press and magazines. 
Mr. Aaron agreed to accept the Chairmanship of this Membership 
Subcommittee on Publicity. 


All Chapter Chairmen are hereby requested to send informa- 
tion that can be used for publicity to Mr. Aaron. A copy of 
papers of interest to the profession in particular and to the 
public in general, should be forwarded in advance of their pre- 
sentation with their date of presentation. This will enable 
releases to technical magazines and the Press. Releases after 
the presentation are almost never accepted. 


Papers Procurement and Publications Committee: It was proposed that 
itor Axelby c € activities of a combined papers study, 
papers procurement, papers review and publications committee. 
Mr. Axelby was authorized to appoint such other members to this Com- 
mittee as may be necessary to carry on its activities. This should 
include the appointment of one person who. would be responsible for 
coordination with other Professional Groups. It is hoped that the 
Abstracts and Reviews Section of the Transactions can be handled on 
an international basis. 


It was the concensus of opinion that the Group sponsor an annual sym- 
posium along the lines suggested by Mr. Axelby as soon as this is 
practicable. After a lengthy discussion it was decided to table this 
point until the next meeting at which time a Symposium Committee will 
be appointed. 


4. Chapter Reports: 


ha. 
lib. 


hee 


hd e 


Los Angeles Chapter: Mr. Grabbe reported that the Los Angeles Chapter 
had suscocded In Increasing its meeting attendance by presenting 
papers on new advances. 

Twin Cities Chapter: Chapter Chairman La Hue reported that the Twin 
Cities Chapter fad enjoyed very good attendance at its meeting on 
"Inertial Navigation." 


Akron Cha ter: It was moved that the petition presented to Headquar- 
ters for tne formation of an Akron Chapter be approved (unanimous). 


New York Long Island and Northern New Jersey: An effort will be made 
Eo contact organizers in these Sections in order to form Chapters. 
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he. Baltimore Chapter: Mr. Axelby reported that two meetings have been 
planned for The coming year and the papers presented may be published. 
in the Transactions. The AIEE Baltimore Section has approached the 
PGAC Baltimore Chapter with the suggestion that joint meetings be 
held. 


Membership Survey: It was decided to prepare and distribute a questionnaire 
to the Membership to determine which of the members are interested in doing 
Committee work, either on a local Chapter level or on a National level, or 
for symposia or papers procurement and review. This questionnaire should 
also attempt to determine what type of programs, meetings and papers for pub- 
lication would serve the membership in the best manner. ‘Fhe Chairman and 
Secretary will prepare the text for this questionnaire. It was suggested 
that in some cases it would be desirable to address letters to Management to 
request their assistance in allowing time for engineers to work on PGAC pro- 
jects. This approach will be followed after the survey has been completed. 


Finance Report: The balance in the Group Tréasury on August 31, 1956, was 


9 e e 


Membership: As of August 31, 1956, membership was as follows: 


Paid 1,482 
Paid Students "157 
Unpaid 159 
Total 1,798 


Date of Next Meeting: The next meeting will be held on Monday, December 35 
1956, at 10:00 a.m. at IRE Headquarters, New York. 


There being no further business, the meeting adjourned at 3:30 p.m. 


M. Robert Aaron 
Secretary-Treasurer 


Issued: September 28, 1956 
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Gitten, L. J. 
Goodstein, Julian 
Gretz, R. W. 
Gross, S. H. 
Hansen, H. R. 
Harris, H. I. 
Haynes, N. M. 
Heacock, W. Je, Jr. 
Heimer, R. L. 
Herman, Sidney 
Hittner, Charles 


. Hoffer, L. H. 


Jacob, G. W. 
Joline, E. S. 
Kay, L. Me 

Kiel, die He 
Klein, R. C. 
Knocklein, H. P. 
Knox, R. W. 
Kokkinos, Constantinos 
Lemanczyk, J. C. 
Levinson, Emanuel. 
Lutz, C. H. 
Marston, R. S. 
Match, M. J. 
McPherson, D. L., III 
Meirowitz, R. L. 
Mohr, H. F. 
Moritz, F. G. 
Morse, R. V. 
Odessey, P. H. 
Orford, R. J. 
Osder, S. S. 
Padovano, R. Jd. 
Pagano, C. N. 
Pearsall, C. H., Jr. 
Peterson, H. 0. 
Pighi,. L. H. 
Prichard, J. S. 
Purdy, Ris. ine 
Rehberg, C. F. 
Rosenthal, M. M. 
Rosenthal, S. A, 
Rubenson, J. G. 
Sandock, Jules 
Sayer, George 
Scheideler, G. E. 
Schneider, R. M. 
Schoedel. R. A. 
Schramm, M. Wes Jr. 
Schulkind, Donald 
Scott, J. E. 
Seckler, P. J. Ac 
Shanahan, W. J. 
Sherry, L. I. 


Sieminski, Edward 
Simon, R. L. 
Simonton, L. J. 
Slusarchyk, G. J. 
Smilowitz, S. N. 
Steinberg, W. A. 
Stephenson, J. G. 
Stepper, L. H. 
Tamasi, A. P. 
Torn, L. J. 
Varmm, A. M. 
Vogel, Erwin 
Walker, A. ©, 
Wathen, R. L. 
Weintraub, Irving 
Wellman, F. C. 
Westover, T. A. 
Wheeler, H. A. 
Wiesner, Leo 
Willis, P. A. 
Winzemer, A, M. 
Wolkov, David 
Wurman, Gustave 
Young, V. J. 


New York Section 
Kibanese, A. P. 
Antonazzi, F. J. 
Arcidiacono, Joseph 
Banner, Leonard 
Barker, D. R. 
Baumann, D. A. 
Beeh, R. C. M. 
Belluck, R. E. 
Bertram, J. E. 
Bibbero, R. J. 
Blecher, Franklin 
Blumberg, R. H. 
Boyanovitch, Dushan 
Boylan, Le Ke 
Brailey, M. L. 
Buton, E. C., Jr. 
Carpenter, E. J., Jr. 
Clemens, G. J. 
Connelly, J. J. 
Cooper, P. C. 
Cristina, A. A, 
Cuccia, J. F. 
Cypser, R. J. 
D'Amato, R, Je 
Dellenbangh, F. S., III 
DeWitt, R. G. 
Dickstein, S. Re 
Diebold, J. T. 
Diehl, Hee To 
Drescher, Charles 
Duffy, J. J. 
Ebaugh, Des Pte 
Evans, B. O. 
Fabricant, B. S. 
Feinerman, Bernard 
Fernald, 0. H. 
Fisch, Alan 

Foran, Re Ee 
Franklin, Ge Fe 
Freudenberg, Boris 
Friedensohn, George 
Gardner, R. Ke 
Garman, R. L. 


Gayet, P. J. 
Georgallis, G. C. 
Gilman, G. We 
Gister, Stanley 
Glazer, Eli 
Glowalla, John 
Goldberg, Harold 
Goldsmith, A. N. 
Goldstein, Edward 
Gronner, Alfred 
Groudan, A, I. 
Haddad, R. A. 
Hamer, Howard 
Heineck, A. We 
Hoekstra, Robert, Jr. 
Hornung, S. A. 
Horowitz, I. M. 
Humphrey, R. M. 
Jorysz, Alfred 
Kalogeras, G. Ae 
Kaplan, K. R. 
Karmel, P. R. 
Kassel, Aaron 
Katz, M. D. 
Klein, W. T. 
Kleinrock, Leonard 
Korrol, ©. R. 
Kram, Bernard 
Kurzrok, R. M. 
Kuzmyak, M. G. 
Lampert, Leon 
Lannary, John 
Liccardi, A. L. 
Lindner, N. J. 
London, F. H. 
Low, Frank 
Maedel, G. F. 
Maitra, K, K. 
Malina, Meyer 
Marcinkowski, H. L. 
Marolda, E. A. 
Marshall, S. L. 
Meyers, S. J. 
Milenbach, Hermane 
Molnar, R. J. 
Ortmann, M. W. 


Otis, A. N., Jr. 


Paschetto, E. Jd. 
Piccione, N. A. 
Porter, R. W. 
Preiss, R, de 
Rabinowich, Josef 
Ragazzini, J. R. 
Redmond, K. P. 
Reeves, J. F. 
Roberts, R. P. 
Rosaler, R. C. 
Rosenberg, A. E. 
Sabo, J. D. 
SanGiovanni, Carlo, Jr. 
Sarachik, Po ne 
Savage, J. E. 
Sherman, Seymour 
Shinners, S. M. 
Slavin, M. J. 
Smith, M. V. 
Snyder, J. B. 
Stanton, L. J. 
Stone, E. 5S. 
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Svigals, Jerome 
Thomas, R. 0. 
Trachtenberg, Alfred 
Aparinells GhG (Cr 
Turchiano, M. We 
Turczyn, W. A. 
Wallace, V. L. 
Walton, J. S. V. 
Warner, F. N. 
Watkins, J. E. 
Watson, E. F. 
Weitman, Irwin 
Wernick, J. I. 
Wolfson, Richard 
Wood, A. G. 


Northern New Jersey Section 
Aaron, M. R. 


Acker, J. L. 
Anderson, WN. E. 
Bahls, W. E. 
Bearman, A. L. 
Bogner, Irving 
Brown, A. T., III 
Brown, C. S. 
Brown, R. I. 
Bucher, F. X. 
Bunko, Myron 
Carr, G. W. 
Cowles, W. W. 
Davis, E. S. 
Doba, Stephen, Jr. 
Doniger, Jerry 
Dorros, Irwin 
Farber, D. J. 
Garfinkel, B. D. 
Glaser, W. A. 
Grandmont, P. E. 
Greenberg, Howard 
Hamming, R. W. 
Helgeson, B. P. 
Huang, R. Y. 
Hunter, W. B. 
Kaiser, E. T. 
Kelly, R. J. 
Klarman, K. J. 
Kulik, V. A. 
Kundel, E. A., Jr. 
Lazos, N. J. 
Loumeau, R. F. 
Lozier, J. C. 
Lunney, R. E. 
McCrory, J. R. 
Moll, John 
Mount, Ellis 
Mueller, P. Le 
Panter, P. F. 
Podell, R. L. 
Reilly, R. A. 
Richards, G. P. 
Ronne, Jo Se 
Russell, F. A. 
Seckler, H. N. 
Seeb, Stanley 
Seergy, C. M. 
Shapiro, Oscar 
Sherin, David 
Sippach, ne Wes Jr. 
Smith, E. J. 


Spock, R. E. 
Streeter, T. W., Jr. 
Sutton, R. H. 
Sweeney, W. R. 
Thompson, C. F. 
Turkheimer, P. M. 
Wall, G. B., III 
Warden, F. W. 
Wilde, A. E., Jr. 
Yamagami, Yonehisa 
Zayac, F, R, 
Zimmerman, Leonard 


Princeton Section 
Hellstrom, M. J. 
Hoedemaker, R. W. 
Lemelson, J. H. 
Losher, M. I. 
Rozanski, R. R. A. 
Ruble, G. B. 
Schofield, C. R. 
Selinsky, J. J. 
Surber, We H., Jr. 
Truitt, T. D. 


REGION III 


Atlanta Section 
Coleman, B. K. 
Eckel, J. Re, Jr. 
Glaser, H. I. 
Lowman, B. L. 
Pippin, R. F., Jr. 
Williams, B. ©, 
Ziegler, N. F. 


Baltimore Section 
Ausfresser, H. D. 
Aviles, H. N. 
Axelby, G. S. 
Bastow, J. G., Jr. 
Behm, G. T. 
Black, F. R. 
Brodwin, M. E. 
Buchan, J. F. 
Choksy, N. H. 
Cichanowicz, H. J. 
Clarke, D. R. 
Coppel, J. M. 
Coulter, E. L. 
Dietz, J. H. 
Edwards, R. L., Jr. 
Fegely, W. D. 
Fuchs, A. M. 
Gambrill, R. D. 
Glaser, E, M. 
Gray, A. R. 
Groszer, A. J., Jr, 
Haut, 0. G., 201i 
Horn, R. E. 
Hurley, W. A. 
Ichniowski, F. ©. J. 
Jentilet, Adam 
Jones, IL. G F. 
Jones, W. N. 
Kegel, A. G, 
Kernan, Paul 
Kintner, P. M. 
Leahy, F. N. 


Muller, J. F. 
Myers, F. G. 
Patton, H. W. 
Penabad, Joseph, Jr. 
Poston, M. H. 
Randolph, G. W. 
Raynes, H. D. 
Roediger, F. E, 
Sovill, A. W. 
Spink, P. G. 
Stebbins, W. J. 
Stefan, R. S. 
Taragin, Saul 
Thomas, Ae ale 
Titen, Harvey 
Tomlinson, C. C. 
Wascavage, J. A. 
Watts, H. M. 
Wilson, H. C. 
Wolf, H. S. 
Wolpert, M. L. 
Zacharia, Harry 


Central Florida Section 
Dibble, H. L. 

Grant, W. S., Jr. 
Howard, T. B. 

Jahimiak, Roger 


Huntsville Section 
Bradley, B. C. 
Cutting, Elliott 
Hatcher, C. R. 
Hickey, R. E. 
Ingram, J. B. 
Schwab, W. G. 

Tu, J.C. 


Miami Section 
Orihuela, Adolfo 
Payne, V. E. 
Rosenzvaig, J. E. 


North Carolina-Virginia Section 
Andrews, R. EB. . 


Bradfute, G. A., Jr. 
Cockrell, W. D. 
Cooper, Benjamin 
Eller, J. E., Jr. 
Gregory, C. A., Jr. 
Imus, H. 0., Jr. 
Jones, T. E. 
Lindeman, W. D. 
Passera, A. L. 
Stephens, T. L. 
Young, D. B. 


Northwest Florida Section 
ollins, J. 0. 
Gamel 9 We. We 


Philadelphia Section 
Affel, H. A., Jr. 


Aires, R. H. 
Anderson, W. G. 
Bachofer, H. L. 
Beaumariage, D, C. 
Beck, Cyrus 
Benner, R. H., II 
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Berg, N. E. 
Booth, A. T., Jr. 
Boyd, W. L. 
Boys, H. N. 
Brandt, W. E. 
Brucklacher, J. E., Jr. 
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Caplan, D. I. 
Carpenter, R. A. 
Cecala, J. A. 
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Cohen, B. H. 
Cole, Je C., JY. 
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Gregory, T. R. 
Hom, F. M. 
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Lee, F. F. 
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Smith, D. B. 
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Sorkin, C. S. 
Stephenson, J. M. 
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Sun, Hun-Hsuan 
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Turner, L..C. 
Tweet, B. 0., Jr. 
walker, H. R. 
Weiner, J. R. 
Weisenberger, A. J. 
Wills, W. P. 
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Wolin, Louis 
Wolin, Sarmel 
Yamada, Hisao 


Washington Secton 
Kiften, Deke 
Britton, D. D. 
Bryant, F. B. 
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Bush, G. B. 
Caggiano, V. J. 


Calhoon, T. G. Furukawa, M. M. Ho, Yu-Chi 


Carruth, D. E. Georger, L. J. S 
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Lee, A. M. 
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Poland, W. B., Jr. 
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Rogers, A. L. 
Rosenzweig, M. S. 
Sanborn, G. D. 
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Singer, J. R. 
Stapleton, J. F. 
Stoops, C. W. 
Varela, A. A. 
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White, C. F. 
Wimmer, P. L. 
Young, H. D. 
Zastrow, K. D. 


REGION IV 


Akron Section 
Colletti, Nicasio 
Diamantides, N. D. 
Flowers, H. L. 
Haas, D. L. 

Hann, D. D. 
Lambert, C. 0. 
Michaelis, T. D. 
Miller, J. H. 
Murray, P. We 
Penniman, tebe 
Rose, We A. 
Ryburn, ye W. 
Sabol, R. W. 
Stahl, M. D. 
Toman, W. J. V. 
Yochelson, S. B. 


Cincinnati Section 
Baird, ° W, 

Berg, De Le 
Colclaser, R. A. 
Dale, W. L. 
‘Doerr, W. H. 
Engelmann, R. H. - 


Willsey, R. H 


Winkeljohann, Albert 


Cleveland Section 
Auth, L. V., Jr. 
Craacrat on Lis 
Dambach, R. A. 
Frenk, D. R. 
Gogia, J. K. 


Grasson, Walter, Jr. 


Hart, C. E. 
Hotchkin, E. E. 
Huntley, J. R. 
Kinkaid, J. C. 
Klock, H. F. 
Mergler, H. W. 
Praff, bs We 


Phillips, W. E., Jr. 


Post, R. H. 
Saltzer, Charles 
Seaton, Glem 
Tame, J. S. 
Trinkle, F. J. 


Columbus Section 
Burgener, R. U. 
Chope, H. R. 
Cohen, Donald 
Conlon, R. J. 
Cook, E. E, 
Fenwick, We D. 
McFarland, R. S. 
Spergel, Philip 
Weimer, F. C. 
Wilson, W. A. 


Dayton Section 
Bornhorst, K. F. 


Charbonneaux, W. A. 
Finnigan, R. E. 
Herrin, C. B. 
Kiebert, M. V., Jr. 
Lewis, D. E. 
Martino, J. P. 
Pastori, D. F. 
Simopoulos, N. T. 
Spengler, J. R. 
VarWechel, R. J. 
Wichmam, T. F. 


Detroit Section 
Barcus, Ronald 
Brown, L. R. 
Bublitz, A. T. 
Chow, Henry 
DaRoza, F., G. 
El-Melehy, M. A. 
Gaskill, R. A. 
Genzlinger, Vance 
Gilbert, E. 0. 
Gilbert, E. G. 
Halsted, L. R. 
Heckler, R. R. 
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Nixon, J. D. 


Olson, RK. G. 
Rauch, L. L. 


Sattinger, I. J. 
Scott, D. E. 
Seleno, A. A. 
Sims, R. C. 
Smith, Wray 
Sutton, W. A. 
Taplin, L. B. 
Theodoroff, T. J. 
Tubbs, R. J. 


Webber, R. C. 
Wyble, J. J. 


Emporium Section 
Bennett, P. E, 


Haines, H. M. 
Harvey, H. B. 
Hoechner, I. L. 
Lemley, L. W. 
Oblinger, J. T. 
Pryima, R. M. 
Seeley, R. M., dr. 


Pittsburgh Section 
Barnard, J. D. 
Brendle, T. A. 
Caywood, W. P., Jr. 
Ellison, B. P. 
Fulmer, L. C. 
Gocsik, J. M. 
Halprin, L. H. 
Kinder, H. R. 
Knowles, C. R. 
O'Donnell, J. J. 


Rau, F. Ale 
Rogers, L. J. 
Royer, G. H. 


Starbuck, W. H. 
Strull, Gene 
Sze, T. W. 


Werst, M. C. 


Toledo Section 
Ewing, D. J., Jr. 
Spademan, C. F. 


Williamsport Section 


Webb, H. Ee 
REGION V 
Cedar Rapids Section 
Hedgecock, W. Te, Jr. 
Lowenberg, E. C. 
Chicago Section 
Antonelli, D. R. 
Axelrod, L. R. 
Bergen, H. A. 


Bold, N. T. 
Buchta, J. C. 


Bullen, C. V. 
Bymberg, R. Je, JYe 
Carter, Robert ~ 
Cermak, C. W. 
Chang, Bansun 
Chulsky, Isadore 
Cooney, J. Jd. 
Dolce, S. L. 

Druz, We os 
Dunbar, E. A. 
Epley, D. L. 
Ferre, G. Ee 

Fu, King-Sun 
Gerlach, A. B. 
Greenberg, C. J. 
Gregory, E. C. 
Hansen, A. G., Jr. 
Hoffman, C. H. 
Hori, Shizuo 
TIsolampi, G. E. 
Jenness, R. R, 
Kott, We O- 

Kreer, J. B. 

Kuhn, N. J. 
Lafferty, V. C. 
Lee, D. K. K. 
Leth, T. R. 

Lewis, H. A. 

Li, Ching-Chung 
Martin, J. W., Jr. 
Merz, R. A. 

Meyer, Andrew 
Mitchell, F. R. 
Mittelmamn, Eugene 
Nauer, M. J. 
Noble, D. S. 
O'Neill, R. M. 
Saltzberg, Theodore 
Shewan, William 
Shively, R. R. 
Smithana, D. R. 
Stan, John 
Thielen, L. R. 
Vaicunas, A. A. 
VanBosse, J. G. 
VanNess, J. E. 
VanValkenburg, M. E, 
Verbanec, W. Re 
Warshawsky, Jay 
Weissert, R. K. 


Evansville-Owensboro Section 


Kercher, D. L. 


Fort Wayne Section 
Ackworth, D. L. 
Brady, F. H. 
Emery, R. C. 
Kalish, J. H. 
Kaplan, Robert 
Norris, B. J. 
Richeson, W. E,. 
Solomon, R. M. 


Indianapolis Section 
Billheimer, A. C. 


Evans, Bly Ae 
Garofalo, A. D. 
Grant, M. P. 


Hammond, S. B. 
Longren, W. K. 
McCrocklin, R. E. 
Ogborn, live slits 
Pearson, Se J. 
Simian, Wo Go the 


Louisville Section 
Kwo, T2 Hib 
Means, T. S., Jr. 


Milwaukee Section 
Arakelian, G. P. 
Babladelis, George 
Cork, H. A. 
Gessner, Urs 
Goodman, P. H. 
Graham, J. D. 
Haraldsen, H. P. 
Jensen, K. 5. 
Lind, E. R. 
Makela, L. V. 
Mezger, Je P. 

Min, H. S. 

Morin, D. C., Jr. 
Pierce, R. L. 
Rekoff, M. G., Jr. 
Schlager, K. J. 
Schwartz, E. B. 
Smith, C..C. 
Zelazo, N. K. 


Twin Cities Section 
Adams, G. i, 
Adkisson, W. M. 
Alderson, R. C. 
Alfsen, G. F. 
Allen, D. H. 
Anderson, L. T. 
Balzart, E. Je, Jr. 
Bartlett, V. We 
Beaudoin, P. E. 
Benassi, D. Ae 
Bergan, K. N. 
Bock, E. D. 
Brunetti, Cledo 
Carlson, R. A. 
Clark, R. N. 
Croze, M. We 
Cummings, K. C. 
Dundovic, J. F. 
Dunwell, R. D. 
Fox, A. J. 
Gilson, J. R. 
Gise,o Pw G. dr. 
Gustafson, H. A. 
Hardenbergh, G, A. 
Hid atom 
Hulstrand, B. E. 
Inman, T. F. 
Johnson, H. W. 
Kershaw, J. A. 
Ketchum, J. R. 
Koenig, J. D. 
Lahue, P. M. 
Lanzkron, R,. W. 
Lindemann, A. W. 
Lode, Tenny 
Ludwig, J. T. 
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Macomber, G. R. 
Markusen, D. L. 
Maze, Re. O. 
McLane, R. C. 
Moe, W. J. 
Muckenhirn, O. W. 
Murphy, G. J. 
Nellis, W. M. 
Peatman, J. B. 
Rowland, C. A., Jr. 
Schuck, O. H. 
Sear, A. W. 
Senstad, P. D. 
Stewart, W. #. 
Stone, N. T. 
Swanlund, G. D. 
Toth, D. H. 
Vogel, J. P. 


Alamogordo-Holloman Section 


Bauman, #. J. 
Kuerschner, Helmut 
Zimmerman, A. P. 


Dallas Section 


—Askew, W. J., JY. 


Braun, Clarence 


Buehrle, C. D. 


Buzard, R. S. 
Creager, L. D. 
Harmon, F. I. 
Heizer, K. W. 
Johnson, G. D. 


-McDonald, Marshall 


Miller, N. D. 
Morris, B. V. 
Ocnaschek, F. J. 
Pitomans be) Des die 
Prier, H. W. 
Stanton, A. N. 
Tatum, F. W. 


. Wadel, L. B, 


Denver Section 
Daniels, W. H. 
Bimen Men Ue 
Messler, F. J. 
Mielziner, Walter 
Morris, Mich Lis 
Morroni, D. J. 
Ostwald, L. T. 
Schneebeck, D. A. 
Tary, J. J. 


El Paso Section 
Emerling, R. A. 
Rojas H, A. M. 


Fort Worth Section 
Beckman, W. R. 
Blasingame, J. J. 
Cone, J. H. 

Evans, W. L. 

Heizer, L. E. 

Jiles, C. W. 

Lowrie, G. M. 
Teasdale, A. R., Jr. 


Watkins, 0. E. 
Young, F. W. 


Houston Section 
Bucy, J. Fey Jr. 
Easterling, M. F. 
Francis, L. G. 
Frobese, C. W. 
Gentry, Je Ae 
Hutchens, R. L. 
Navarro, 5S. 0. 
Tasini, Betsalel 
Waldie, A. D. 


Kansas City Section 
Breyfogle, L. D., III 
Hickey, L. F. 

Miller, H. G. 

Murray, W. A. 

Stout, H. L. 

Wilcox, J. V. 


Little Rock Section 
annon, W. W. 


Lubbock Section 
Estes, 5S. KE. 
Perkins, C, S. 
Tomlinson, Z. G. 


New Orleans Section 
Drake, R. L. 


Oklahoma aie Section 
Grubbs, C. 4. 
Krystek, M. E. 
Ledbetter, R. P. 
Puckett, T. H. 

Silva, R. F. 

Viay, G. J. 


St. Louis Section 
Allison, W. H. 
Arndt, R. L. 
Hibbits, R. M. 
Lago, G. V. 
Malsbary, J. S. 
Mayer, M Es, Jr. 
Mohrman, R. F. 
Mutchek, J. H. 
Norman, C. F. 
Reed, D. L. 

Salman, N. D. 
Sayer, J. D. 
Tedeschi, Anthony 
Twombly, J. W., Jr. 
Weinstock, G. L., Jr. 
Winter, Demle 


San Antonio Section 
ostick, F. X. 
DuBose, G. P., Jr. 
|Hirsch, C20 
‘Hoffman, A. A. J. 
ayleben, E. F. 
Reinhard, E. A. 
Ziemer, D. R. 


Tulsa Section 
Brashear, Rk. T. 


Cairns, T. W. 

Day, C. E, 

Fox, D. N. 
Labarthe, L. C. 
Laird, J. Aw, Dif 
O'Brien, D. G. 
Rowley, R. G., Jr. 
Silverman, Daniel 
Sykora, G. E. 


REGION VII 


Albuquerque-Los Alamos Section 
Ehrman, ocd 


Katzenstein, Jack 
Pace, T. L. 
Shephard, R. W. 


China Lake Section 
Crawford, J. A. 
Creusere, M. C. 
Kim, P. K. S. 
Schimmel, George 


Hawaii Section 
Jones, R. C., Jr. 


Los Angeles Section 
Abbott, W. R. 


Ackerlind, Erik 
Akin, P. A. 
Albrecht, Albert 
Anderson, Frank 
Andrews, L. A. 
Anzel, B. M. 
Aroyan, G. 1 
Aseltine, J. A. 
Avrech, Norman - 
Baker, D, L. 
Ballard, K. C. 
Barlett, F. R. 
Barnes, J. L. 
Bayley, L. B., Jr. 
Beckwith, H. W. 
Beecher, A. E. 
Bekey, G. A. 
Bement, W. A. 
Bible, R. E, 
Bills, G. W. 
Bonney, R. B. 
Borgeson, P. W. 
Bower, J. L. 
Braverman, D. J. 
Broadwell, W. B. 
Brock, P. A. 
Brown, D. E. 
Buchman, W. W. 
Buland, R. N. 
Burk, W. A. 
Burnsweig, Joseph, Jr. 
Carlson, A. R. 
Carlson, C. 0. 
Carlson, R, S. 
Cassidy, R. E. 
Chandaket, Prapat 
Chandler, D. P. 
Christensen, A. V. 
Cordray, R. E. 
Corvi, J. A. 
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Cosgrave, S. J. 
Curry, W. S., Jr. 
Deaux, F. J. 
Deming, A. F. 
Deuser, D. A. 
Dickinson, H. B. 
Diem, C. W. 
Diemer, F. P. 
Dinning, J. R. 
Doty, R. L. 
Drucker, Alvin 
Dzilvelis, A. A. 
Edelsohn, C.'R. 


Eikelman, J. A., Jr. 


Engel, H..L. 
Eno, R. F. 


Fernandey, Ferdinand 


Finley, W. A. 
Fish, W. Y. 
Forbath, F. P. 
Foxman, Eugene 
Francis, T. F. 
Frankel, Sidney 
Frankos, D. T. 
Fuller, R. H, 
Fulton, A. 5. 
Furumoto, Nobuo 
Gabler, R. T. 
Gaitan, T. C. 
Garber, L. F. 
Gauronskas, P. P. 
Gaylord, R. S. 
Gerardi, F. R. 
Gerken, G. H. 
Gill, W. J. 
Ginstling, Ajzyk 
Grabbe, E. M. 
Graham, J. D. 
Gross, William 
Gunning, W. F. 
Hadden, F. A. 
Hansen, 0. B. 
Harmon, W. G. 
Harrington, L. M. 
Hassel, R. R. 
Hayes, J. E. 
Heffner, E. K, 
Heyliger, G. E. 
Hicks, A. R. 
Hruby, R. J. 
Hughes, Frank 
Hutcheon, R. S. 
Izuel, A. G. 
Jack, R. W. 
Jackson, K, R. 
Jacobs, Je E. 
Janeway, R. K. 
Joerger, J. C. 
Johnson, J. Jd. 
Johnson, R. We. 
Johnson, W. A. 
Juran, Warren 
Kaufman , Ke ie 
Kaufman, Sidney 
Kawahata, B. I. 
Kennel, J. M. 
Keppel, R. A. 
Kerster, George 
Waiteigs Gg Cig dhe 


King, J. E. 
Kirsch, H. A. 
WGislasbs 15 Jol 
Klein, M. L. 

Knox, R. V. 
Krames, C. V. 
Geshe, 105 Is 

Kroy, W. H., Jr. 
Lawrence, A. F., III 
Lee, H. Jd. 
Leondes, C, T. 
Leone, W. C. 
Levinson, R. M. 
Lewis, D. E. 
Liang, Ming-Tsu, M. 
Lillibridge, E. H. 
ViOTSis Nie, Vis 
Louie, William 
Lyons, L. H. 
Malone, Martin 
Mancini, A. R. 
Mankinen, Ee J. 
Manly, Ron 
Margolis, Maier 
Mayberry, L. A. 
McRuer, D. T. 
Mehner, E. W. 
Mileson, D. F. 
Miller, D. S. 
Mitsutmoi, Takashi 
Morrison, A. I. 
Morton, W. B., Jr. 
Myers, W. A. 
Nelson, C. S., Jr. 
Neumann, Leopold 
Noland, A. R. 
Nuban, Ebrahim 
Nuttall, H. V. 
O'Brien, W. C. 
Olsen, L. V. 
Parker, A. T. 
Pernick, L. J. 
Post, Geoffrey 
Poulson, W. A. 
Primozich, F. G. 
Putter, Klaus 
Quackenbush, R. E. 
Radant, M. E. 
Raffensperger, M. J. 
Ramer, F. H., Jr. 
Ramstedt, C. F. 
Redden, E. T. 
Redmond, J. G. 
Rehler, K. M. 
Rescoe, J. M. 
Rickords, T. J. 
Rifkind, Jesse 
Robertson, G. R. 
Rogers, J. G. 
Rogers, T. A. 
Romano, A. J. 
Rosenstein, A. B. 
Rosenthal, G. W. 
Rowe, D. E, 
Ruiz. Le 
Salzer, J. M. 
Samuels, A. H. 
Sanneman, R. W. 
Sarture, C. We 


Savant, C. J., Jr. 
Savo, T. A. 
Sawyer, H. F. 
Scammel, B. C. 
Schalk, Norbert 
Schroeder, William 
Schultz, P. R. 
Schultz, R. T. 
Schulz, K. 5. 
Scott, W. F. 
Sensiper, Samuel 
Shelley, R. G. 
Shenk, J. We 
Shimada, George 
Short, F. E. 
Shuler, M. H. 
Shultise, Q. M. 
Shutt, S. G. 
Siegel, J. C. 
Silva, L. M. 

Sink, R. L. 
Slocomb, G. M. 
Smith, J. C. 
Smith, J. E. 
Snapp, K. M. 
Snyder, W. Ae 
Sohler, J. F. 
Stark, H. P. 
Steinkolk, R. B. 
Stimpson, L. D., er. 
Stout, Ue Me 
Sturm, T. F. 
Takahashi, Kiyoshi 
Thompson, D. M. 
Thorensen, Ragnar 
Thorpe, L.—M. 
Turn, Rein 
Valery, N. A. 
VanCuren, Verlyn 
Vega, C. d. 
Vittum, W. M. 
Vodovoz, Erwin 
Vulliet, P. 0. 
Wachowski. H. M. 
Waddell, k. L. 
Wakamiya, Yooichi 
Walker, N. L. 
Walkup, L. A. 
Wallace, Charles, Jr. 
Waln, R, M. 
Walters, L. G. 
Wanlass, S, D. 
Warrington, Wilmer, Jr. 
Watkins, E. L. 
Wedel, J. J., Jr. 
Wennerberg, Gunnar 
Wenters, R. L. 
White, L. M. 
Whitford, R. K. 
Wolman, L. L. 
Wong, D. S. 

Wong, E. C. 

Wong, Herbert 
Young, W. L. 
Zabusky, N. J. 
Zacharias, Robert 
Ziegler, R. M. 
Zimmerman, R. L. 
Zoller, ©. J. 
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Phoenix Section 
Ballantine, J. H. 
Levine, Daniel 
Ross, J. M. 
Scrafford, R. L. 
Sutton, J. G., Jr. 


Portland Section 
Bowman, D. G, 
Doel, Dean 
Jenkins, R. W. 
Mendoza, D. C. 
Stone, L. N. 


Salt Lake City Section 
Clegg, J. C. 


Gilmour, G. A. 
Jackson, H. L. 
Murphy, L. C. 


San Diego Section 
Biering, A. H. 


Campbell, Robert 
Cox,a he Ms 

Dodd, G. M. 
Ferner, R. 0. 
Flarity, W. H. 
Fogel, L. J. 
Greenspan, L. E. 
Kalbfell, D.C. 
Klimberg, Joseph 
Lay, P. D. 
Mealey, G. J. 
Okahata, E. S. 
Shechet, M. L. 
Wade, Ernest 


San Francisco Section 
Kiexander, C. S., Jr. 
Berryhill, J. L. 
Bharucha, B. H. 
Binnall, E. P. 
Brennan, R. D. 
Buntenbach, R. W. 
Chesebro, E, L. 

Chin-e bens 

Davy, L. H. 

Dibb, George 

Durfey, G. K. 

Enos, R. M. 
Farman-Farmaian, Gheffer 
Firschein, Oscar 
Gardiner, K. W. 
Gentry, E. B, 
Gleason, C. A. 
Goodall, J. R., Jr. 
Gyllstrom, N, D. 
Hallman, A. B. 

Hexem, John 

Iwama, Morimi 
Jameson, R. J. 

Jury, Merle 
Kazanjian, H. A. 
Kerwin, W. J. 

Kincas, J. W. 
Klotter, Karl 
Kochenderfer, W. E., Jr. 
Kortman, C. M. 
Landsman, Louis 


Lendaris, G. G. 
Lessley, T. D. 
Lindberg, H. E. 
Martinez, H. M. 
MUU yee Berd. 
Ness, R. J. 
Nishizaki, Ray 
Oliver, BoM: 
Over, J. J., Jr. 
Pope, J. C. 
Pringle, Ralph, Jr. 
Rearick, H. F. 
Regenos, K. M. 
Ross, Albert 
Samario, E. J. 
Smith, D. L. 
Smith, 0. J. M. 
Storke, F. P., Jr. 
Tudor, B. J. 
Tuttle, D. F., Jr. 
Vea,ot. H. 

Wang, P. K. C. 
Watson, J. K. 
Windsor, R, N. 


Seattle Section 
Biggs, Je De, Jr. 
Bireh, J. 5. 
Bishop, D. Jd. 
Galloway, W. C, 
Homitch, J. M. 
Ledray, William 
Miller, J. J., Jr. 
Noland, L. Jd. 
Smith, R. A. 
Stapleton, E. R. 
Vermilion, E. E. 
Whipple, M. M. 


Tucson Section 
Bard, W. 4, 

Crow, RB. 
Lindenberg, E. C. 
Martin, L. C. 
Sakrison, D. J. 
Simpson, E. J. 


REGION VIII 


Bay of Quinte Section 
Brule, R. J. 


Flemons, R. S. 
MacKelvie, J. S. 


Hamilton Section 
arnahan, ©. W. 

Kassner, John 

Rogers, A. E. 


London Section 
etcher, H. R. 
Stroud, E. L. 


Montreal Section 
Bar-Urian, Moshe 
Birman, Gerhard 
Caron, J. Y. 
Cummins, J. A. 
Demers, Pierre 


Germain, L. V. 
Ginsburg, Isaac 
Heckman, G. R. 
Malcolm, F. We 
Prichodjko, Alexander 
Reeves, Rene 
Richard, G. B. 
Rinfret, C. J. 
St. Onge, J. L. 
Semple, E, R, 
Sproul, R. W. 


Northern Alberta Section 


Wood, D. R. 


Ottawa Section 
Beneteau, P. J. 
Newton, K, G. D. 


Toronto Section 
Baldwin, J. H. 
Balmain, K. G. 
Beaudoin, D. P. 
Byers, H. G. 
Carew, S. J. H. 
Carley, R. R. 
Elliott, W. F. 
Hackbusch, R. A. 
Lang, G. R. 

Maw, J. C. 
McCloskey, K. P. 
Newhall, E. E. 
Otsuki, J. 5. 
Penrose, R. M. 
Peprnik, H. 0. 
Stoddart, T. W. H. 
Unger, J. H. W. 
Wall, Ernest 


Vancouver Section 
Bohn, E. V. 
Hewit, Hoo. 
Kersey, L. Re 
McDonald, W. H. 
Moore, A. D. 


OVERSEAS MILITARY 
Cassidy, J. J. 
Dunham, P. N. 


FOREIGN SECTIONS 


Buenos Aires Section 
inasco, 5. F. 


Israel Section 
Shamir, Jedidiah 
Weislitzer, Josef 


Rio de Janeiro Section 
DeMattos, H. v. 


Tokyo Section 
Honda, Tatsuo 
Tbuka, Masaru 
Ishikawa, Takeji 
Twakata, Hideo 
Koichibara, Tadashi 
Konomi, Mitsugu 
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Mita, Shigeru 
Morita, Masasuke 
Nakahara, Fujio 
Nishino, Osamu 
Okada, Minoru 
Owaki, Kenichi 
Tanabe, Yoshitoshi 
Tanaka, Yoneji 
Taniguchi, Fusao 
Yano, Akira 


FOREIGN 


Australia 
Brodribbee hcp. 
Davies, R. J. C. 
Honnor, W. W. 


Brazil 
Barros-Barreto, L. A. G. 
Waeny, J. C. GC. 


Denmark 
Jarlov, Ae L. 


England 

Clark, J. E. 

Cullen, A. L. 
Dawes, E. J. 

Dietiker, Walter 
Fleming-Williams, B. C. 
Harris, K. E. 

Laverick, Elizabeth 


-Parsons, A. N. 


Warr, H. J. J. 


France 

Baron, Jean 
Berline, S. D. 
Blachier, B. L. 
Ferrier, P. A. 
Ghertman, Jean 
Girerd, J. L. M. 
Loeb, J. M. 


Germa: 
Busch, C. W. 


_ Peters, J. F. 


Rohde, Lothar 
Uhrmamn, A. M. 


Holland 
Alma, G. H. P. 
Janssen, J. M. Le 
Tellegen, B. D. H. 


India 
Mirchandani, I. T. 
Mishra, Srinibas 
Mukerji, M. M. 


Italy 
Bacchialoni, F. L. 
DeDominicis, C. M. 
Egidi, Caludio 
Pinolini, F. 
Tchou, Mario 
Verdoni, L. G. 
Vergani, Angelo 


Lebanon é Sweden Switzerland — Tiree 
Hoffman, J. D, Andersson, K, N. Braun, Acris ae er 
Mone ay. Ekelof, Stig . Demieville, Henri >? Se 
| Elfving, A. Le Shah, R. R. — 
deepencar shor) ALM, Fagerlind, S. G. Strohschneider, walter ham 
uarez-Diaz, Jorge Gyllenkrok, Thure-Gabriel Thalmann, Victor 
Norwa Josephson, B, A. S. Weber, G. C. Ce aes 
mnevik, S. Be “s — Lofgren, E. 0. 4. 
wee Roll, Anders Venezuela {ae 
Puerto Rico Romeil, G. D. R. Arreaza, R. G. are 
Maldonado, Bolivar Svaia, C. G. Kinzbruner; Paul <i) .s Sees 


es ). 


vp ete os ay 


